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FORTRAN PROGRAM FOR INDUCTION MOTOR ANALYSIS 
by Gary BoIIenbacher 
Lewis Research Center 


SUMMARY 

A FORTRAN program for the analysis of squirrel- cage induction motors is de- 
scribed. The analysis encompasses calculations of torque-speed characteristics, elec- 
trical characteristics, magnetic flux densities, weight, and various other parameters. 
Detailed instructions for use of the program are given. The analysis equations are 
documented, and the sources of the equations are referenced. The appendixes include 
a FORTRAN symbol list, a complete explanation of input requirements, and a list of 
error messages. 


INTRODUCTION 

A FORTRAN program that has been used at the Lewis Research Center to analyze 
the electromagnetic design of three-phase induction motors is described in this report. 
The analysis equations used throughout the program are those commonly found in the 
literature. It is the purpose of this report to document these equations, to cite the 
source references, to provide instructions for using the program, and to assist in inter- 
preting program output. The report also provides information to facilitate program 
modification. 

Specific information regarding the computer program is given in the appendixes: 
Appendix A gives the input requirements. A typical set of input data and the resultant 
output are shown in appendix B. Program listings, error messages, and FORTRAN 
symbol lists are provided by appendixes C, D, and E, respectively. 

The input data and computer output shown in appendix B are for a 1200-hertz induc- 
tion motor designed for the Brayton program. This motor, which operates with the 
cavity filled with oil, has been built and tested at the Lewis Research Center. Prelim- 
inary comparison of the test results with the computer program output shows the com- 
puter analysis to be very accurate. 


The program described in this report uses only U. S. customary units. The Inter- 
national System of Units (SI) is used in the main text to conform to publishing require- 
ments. For clarity in presenting the program, only U. S. customary units are used in 
the appendixes. 


MOTOR TYPES SUITABLE FOR PROGRAM ANALYSIS 

The program described in this report can be used to analyze only a three-phase in- 
duction motor with a squirrel-cage rotor winding. The analysis is limited to steady- 
state, balanced conditions with the motor operating in the motoring mode. 

The armature is assumed to have a two- layer, Y- connected winding. The winding 
material may be either copper, aluminum, or brass. 

The rotor stack is assumed to be the same axial length as the stator stack. Double 
squirrel-cage windii^s are not permitted. Materials for the rotor winding, as for the 
armature winding, may be either copper, aluminum, or brass. Deep rotor bars are 
permitted, but no correction is made in the analysis for nonuniform current distribution 
in the bars at high percentages of sUp. 

A cross-section of a typical induction motor assumed for this analysis is shown in 
figure 1. Rotor and stator slots may be open or partially closed; they may be rectangu- 
lar, trapezoidal, round, or trapezoidal with a rounded bottom. If the slots are trape- 
zoidal, the tooth width may be constant or the slot sides may diverge at any other speci- 
fied angle. 

The effects of windage on motor performance may be either omitted or included in 
the analysis. If they are to be included, the windage loss for a similar motor must be 
known. The program will then scale this known windage loss with respect to all perti- 
nent parameters in order to approximate the windage loss of the motor under analysis. 
The scaling is unnecessary if windage loss at synchronous speed can be supplied to the 
program directly. 

Thermal analysis is not a part of the program. Hence, winding temperatures must 
be calculated or estimated for input to the program. 


PROGRAM DESCRIPTION 
General Information 

The computer program consists of a main program called INDMTR and five subrou- 
tines called CIRCT, MAGNET, SLOTS, WDGFCT, and CMBNTN. The program also 
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calls a data-plotting subroutine PLOTXY (ref, 1), which is part of the Lewis subroutine 
library. 

The main program INDMTR handles all input and output functions (with the excep- 
tion of some error messages) and most of the calculations. Figure 2 shows a simplified 
flow chart of the main program and its subroutine usage. A brief description of each 
subroutine follows. 

Subroutine CIRCT treats the induction motor as a two-terminal electrical circuit 
consisting of resistors and inductors (fig. 3). Analysis of the electrical circuit is equi- 
valent to analysis of the motor and allows the determination of all motor performance 
parameters of interest. Subroutine CIRCT performs the necessary circuit analysis. 

The numerical values of the resistive and inductive circuit elements are calculated by 
the main program. 

Subroutine MAGNET computes the flux densities in the stator and rotor back iron 
and in the stator and rotor teeth. It also computes the ampere-turn drops throughout the 
magnetic circuit. If any portion of the motor saturates, an indicator is set to alert the 
main program. 

Subroutine SLOTS is called twice by the main program: once for the stator slots 
and once for the rotor slots. Its function is to compute all slot dimensions that are not 
input to the program but that are needed in subsequent calculations. SLOTS also com- 
putes slot areas and the slot permeance ratio. 

Subroutine WDGFCT computes the distribution factor and pitch factor for the stator 
winding. It also checks that the winding specified is physically realizable; if it is not, 
an error message is printed. 

Subroutine CMBNTN checks if the number of rotor slots, the number of stator slots, 
and the number of poles are mutually compatible. An incompatible combination is one 
that may result in noise or vibration problems or one that may cause imdesirable torque- 
speed characteristics. If compatibility exists, control is returned to the main program; 
if not, the subroutine prints out an error message and lists the number of rotor slots 
that may be used instead. 

The program is written in the FORTRAN IV programming language for use on the 
Lewis 7044-7094 direct-couple system. On this system, typical preexecution time, in- 
cluding compiling, is 80 seconds. Approximate execution time is 4. 0 seconds for each 
motor analysis (more than one analysis may be performed with each computer run). 

The core storage requirement is approximately 13 000 (decimal) words. 


Synopsis of Input Requirements 

Program input consists of one or more data sets. Each data set contains three 
t 3 T)es of data: 
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(1) Magnetic materials data: These data consist of two material decks, one for the 
rotor and one for the stator. Each material deck contains magnetization curve data and 
core-loss data. The core-loss data may be omitted for the rotor material. 

(2) Windage loss data: These data include a known windage loss at a known (refer- 
ence) condition. The program scales the given windage loss from the reference condi- 
tion to the motor design. A windage loss data card must be included with the input data, 
although windage data may be omitted. 

(3) Motor design data: These data include all physical dimensions that affect the 
motor electromagnetic performance. Also included are the winding temperatures and 
the electrical rating, such as voltage and frequency. Any number of motor design data 
sets may follow the materials and windage data. 

Material data decks must appear in pairs in the input stream. They must be fol- 
lowed by the windage data. Following the windage data are the motor design data decks. 
The material data and the windage data apply to all motor design data decks that follow 
until new material and windage data are encountered by the program. The appearance 
of a new magnetic materials data deck signals the beginning of a new data set. 

To keep ke 3 q)unching to a minimum, much of the input is optional and need not be 
read in. If optional data are omitted from the input, assumptions regarding the omitted 
data are made internal to the program. 

Detailed discussion of the input requirements is included in appendix A. Appendix A 
also identifies the optional data and explains the assumptions made for omitted data. 


Output 

For each motor design data deck the program produces six pages of printed output. 
The first three pages provide a record of most input plus other dimensions or param- 
eters calculated by the program, such as reactances and weights. The next two pages 
give the motor performance for slip values from 1 percent to 100 percent in 1-percent 
increments. Motor performance consists of calculated values of torque, input power, 
input current, power factor, efficiency, and a loss breakdown. The sixth page gives a 
plot of torque against slip. 

Each time a new data set is encountered, one page of output is printed that summar- 
izes the material properties appl 3 dng to all subsequent motor designs. This summary 
consists of magnetization curve data for both the stator and rotor material and core- loss 
data for the stator material only. Appendix B shows a t 3 q)ical output. 

In addition to the standard output, the program prints error messages as necessary. 
Appendix E lists all possible error messages, identifies the subroutine in which each 
message originates, and gives the probable cause of the error. 
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Calculational Methods 


General . - The key to the induction motor analysis is the equivalent circuit shown 
in figure 3. In the circuit, R1 represents the phase resistance of the armature winding, 
XI the armature leakage reactance per phase, and R2 the resistance of the rotor winding 
referred to the armature windii^. The resistance element R2/(S/100) shown in figure 3 
represents the combined effect of the rotor winding resistance and the shaft load. The 
symbol X2 is the rotor winding leakage reactance per phase. Both R2 and X2 are re- 
ferred to the armature winding. Resistance RO and reactance XO allow for the effects of 
core loss and magnetizing current, respectively. Values of Rl, R2, XI, and X2 are 
calculated from the physical dimension of the motor. Values for XO and RO are arrived 
at through an iteration process in the no-load magnetic calculations. The symbol S 
denotes the rotor slip in percent. 

The values of the circuit elements, once calculated, are assumed to be constant for 
all values of slip. Motor performance is computed by doing steady-state alternating- 
current circuit analysis for each pertinent value of slip. 

All equations given are expressed in FORTRAN notation including FORTRAN sym- 
bols, FORTRAN arithmetic operations, and FORTRAN function names. The only de- 
parture from FORTRAN notation is for fractions, which are written upright for better 
readability. 

The SI units are used throughout the main body of the text, hi particular the follow- 
ing units are employed consistently: 

(1) Length in meters 

(2) Areas in square meters 

(3) Resistances and reactances in ohms 

(4) Currents and magnetomotive force in amperes 

(5) Voltages in volts 

(6) Power in watts 

(7) Frequency in hertz 

(8) Rotational speed in revolutions per minute 

(9) Resistivity in ju- ohm- meter 

(10) Mass in kilograms 

(11) Magnetic flux density in webers 

Viscosity may be expressed in any unit provided consistency is maintained. 

Slot permeance ratios . - The rotor and stator slot permeance ratios are required 
to calculate the primary and the secondary slot- reactances. Both slot permeance ratios 
are calculated in subroutine SLOTS. The equations used are discussed in this section. 
Six slot shapes are allowed by the program; the permeance ratio calculations are simi- 
lar for all. 
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Equation (7. 9) of reference 2 gives the slot leakage permeance for an open slot, 
such as slot type 1 (fig. 4) in the program. The equation is repeated here. 




- f) 


( 1 ) 


This equation is rewritten below using the FORTRAN symbols of subroutine SLOTS. 
These symbols are defined in appendix D and in figure 4. 


AXX = * ( D2X + _ D5X * (kx - 0. 6667) (2) 

WSX V 3. / 12.*WSX 4.*WSX 


Rearranging the equation gives 


AXX = KX* 



(3) 


where 


A1 = 0.25*KX + — 

12 . 


and 


A2 = 0.25*KX - — 

6 . 


The ratio D2X/WSX in equation (3) gives the contribution of the slot opening to the total 
slot permeance ratio for an open slot. Thus, the equation may be rewritten as 


AXX = KX* 


/Permeance ratioN 




for slot 
opening 



(4) 


This equation may be generalized to other slot types by substituting the correct ex- 
pression for the slot-opening permeance ratio and by replacing WSX by the slot width 
appropriate for the slot type. The results of these substitutions for the various slot 
types are as follows; 
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For slot type 1, 


axx = kx* Vai-^2L'|.a2 

\wsxJ \wsxJ Vwsx/ 

For slot type 2, 

= 2.*KX*D2X ^ /.D13LUa1 - ( 2.*D5X \ 

WSX + WSX4 \WSX4/ \WSX4 + WSX5/ 


( 5 ) 


(6) 


For slot t3^)e 3, 


AXX = KX* 


D3X 


D4X 


wsxi Vwsx - WSXly 


> . alog * (m .) . A1 - (^'1 ♦ A2 

J Vwsxi/ wsxj Vwsx/ Vwsx/ 


For slot tyi)e 4, 


( 7 ) 


AXX = KX* 


D3X 


D4X 


WSXI VWSX2 - WSXI, 


*ALOGf'^U 


2.*D2X 


WSXI/ WSX2 + WSX4. 


.WSX4/ 




2.*D5X 


WSX4 + WSX5 


*A2 (8) 


For slot type 6, 


AXX = KX* 


D4X 


WSXI VWSX2 - WSXI/ 


Valog 


2.*D2X 


WSX2\ 


wsxi/ WSX2 + WSX4. 


VWSX4/ V 


2.*D5X 


WSX4 + WSX5/ 


*A2 (9) 


For slot type 5 (round slot) the expression for AXX is not readily derived from equa- 
tion (4). Instead AXX is computed as shown on page 235 of reference 2. The equation, 
in FORTRAN symbols, is 


AXX = KX* ( + 0. 625* KX 

\WSXlJ 


( 10 ) 
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The equations just given apply to both the stator and rotor slots. However, for 
rotor slots the factor KX (pp. 184 and 185, ref. 2) is always equal to unity, and the slot 
dimension D5X is assumed to be zero. To change these equations to the FORTRAN 
symbols used in the main program, change the X in each symbol to S for stator slots 
or to R for rotor slots. The only exception is AXX, which changes to AXS and AXR, 
respectively. 

Equivalent circuit elements. - The equivalent circuit elements calculated directly 
from the physical dimensions are Rl, R2, XI, X2, and XOAG. The calculation of cir- 
cxiit elements XO and RO is described in the section No-load magnetic c alculation. 

That component of XO that is attributable to the airgap of the motor is called XOAG. 

It is used as an initial estimate of the value of XO in the iteration procedure employed to 
compute XO. It is also required to calculate several other reactances, notably the skew 
reactance, the rotor and stator zigzag reactances, and the peripheral airgap reactance. 
Since XOAG is a function of the physical dimensions of the motor only, its calculation is 
explained here. 

The quantity XOAG is computed as shown in equation (7. 1) of reference 2. The equa- 
tion is repeated here. 


6. 38 qfN^K^K^DL 

Xj^ = ohms /phase (11) 

X 10® 


To obtain the equation used in the program, the number of phases is set to three and the 
factor kj is set to unity. (In the reference, P is the number of pole pairs and N is 
the number of stator turns in series per phase. In the program, P refers to the number 
of poles, and N is the number of conductors in series per phase. ) Thus, in FORTRAN 
notation and with the sjTnbols of the main program. 


XOAG = 


3.02E-5* ( 

V2.*KPS*KDS/ 


* * 2* D* L* F 


P* P* GE 


( 12 ) 


The armature leakage react ance XI is made up of several components as shown in 
the following equation: 


XI = XSS + XSE + XSK + XSZ + XP 


(13) 
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where 


XSS primary slot leakage reactance 
XSE stator end- connection leakage reactance 
XSK one-haU of the skew reactance 

XSZ stator zigzag leakage reactance 

XP peripheral airgap leakage reactance 

The individual components of the armature leakage reactance are computed as follows: 


XSS = 2. 36E-5* N* N* F* L* 

Vqs/ 

(14) 

9. 45E-6* N* N* F* (KPS* KDS)** 2 * /g , FI ^ DSS\ 
P V 2. 4. / 

(15) 

XSK - 0 5* * ^P* :(e 2 

12. V* D / 

(16) 



XSZ = 0. 833* XOAG* 

(KPS* KDS)* * 2 5 * * * 2 

(17) 

XP = 0. 525* XOAG* *^^ * * 2 

(18) 


For equations (14), (17), and (18), see reference 2 (eqs. (7.2), (7.47), and (7.69), re- 
spectively). For equation (17) also refer to reference 2 (top of p. 202). Equations (14) 
and (15) are documented in reference 2 (pp. 336 and 337). The factor KS is discussed 
on pages 184 and 185 of reference 2. 

The rotor winding leakage reactance X2 is similarly computed. The equations 
used in the program are as follows: 

X2 = XRS + XRE + XSK + XRZ (19) 
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where 


XRS secondary slot leakage reactance 
XRE rotor end- connection leakage reactance 
XSK one-half of the skew reactance 

XRZ rotor zigzag leakage reactance 


XRS = 2. 36E-5* N* N* F* L* 


/AXR\ 
V NB/ 


* (KI^*KDS)**2 


( 20 ) 


XRE = 28. 54* AY * 
P 


2.*P*BR + 


3. 1416* D* DC 


1. 7* TER + 0. 6* (DERI - DER2) + 1. 4* DC. 


( 21 ) 


where 


_ N*N* F*(KPS*KD S)**2*2.4E-7 
P 

XSK = 0. 5* SKEWj * * 2 (22) 

XRZ = 0.833*X0AG* (23) 

(KPS*KDS)* * 2 5 * * * 2 


Equations (21) and (22) are given in reference 2 (p. 237). Equations (20) and (23) are 
described in reference 2 (eqs. (7. 11) and (7.47), respectively). Also refer to refer- 
ence 2 (top of p. 220) for further modification of equation (7. 47). 

The armature re sistance R1 is given by 

_ (LS*N*RSTVT Y* 1. OE-6) (24 

(PC*SS) 

where RSTVTY is the resistivity of the stator winding material at 20° C. Then R1 is 
corrected for the windii^ temperature specified in the input data. 

The rotor resist ance R2 (referred to the armature winding) is the sum of the end- 
ring resistance and the resistance of the rotor bars. It is computed by (eq. (206), 
ref. 3) 
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R2 = [3. * RSTVTY* (N* KPS* KDS)* * 2] * /LB - TER ^ ^64;;^ PgRl* KRING\ ^g5) 

\ SB*NB P*P*SER / 

The end-ring thickness TER is subtracted from the rotor bar length LB because it is 
assumed that the axial c\u*rent flow does not extend to the ends of the rotor bars. The 
factor KRING is included to allow for unequal current distribution in the rotor and 
rings. This factor is fully described in reference 4. Like Rl, the value of R2 is 
adjusted for the specified winding temperature. 

Core-loss calculations . - For all core- loss calculations, the value of WFE must 
be known. This is the core- loss expressed in watts per unit of mass at a given lamina- 
tion thickness, frequency F, and flux density BK. Its value is computed from the input 
data contained on the $FELOSS data cards that are part of the stator material deck. Up 
to ten $FELOSS data cards are allowed with each material deck, one for each lamina- 
tion thickness for which calculations are anticipated. The program searches through the 
FELOSS cards to find the data for the lamination thickness nearest that specified in the 
motor design deck. The core loss WFE is then calculated by 

WFE = WCORE* / ^ ** SLOPE (26) 

Vfcore/ 

where F is the motor design frequency as specified in the motor design deck and 
FCORE, WCORE, SLOPE are given on the FELOSS data cards. (The symbols FCORE, 
WCORE, and SLOPE are defined in appendix A, in fig. 11, and in appendix E . ) Com- 
putation of no-load core loss is deferred until the no-load magnetic calculations. The 
no-load core loss is then used to compute the value of RO in the equivalent circuit. 

Core loss at all other loads is the power dissipated in the resistance RO as determined 
during the equivalent circuit analysis. 

No-l oad magnetic calculations . - The magnetic calculations accomplish several 
things: 

(1) Computation of flux densities throughout the magnetic circuit of the motor at 

no load 

(2) Computation of ampere-turn drops across various parts of the magnetic circuit 

(3) Computation of magnetizing current and magnetizing reactance XO 

(4) Computation of values of core loss WO and the value of the resistive element 

RO used in the equivalent circuit to represent the core loss 
The calculations are performed as follows: Initial extimates for XO, WO, and RO are 
made by using the equations 


XO = 0. 5*X0AG 


(27) 
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wo = 3. * (WSYOKE + WSTOTH)* WFE 


(28) 


R0 = 5-*V1*V1 (29) 

WO 

Next subroutine CIRCT is called to compute the airgap voltage V2 and the magnetizing 
current IMAG (fig. 5). From V2 the total flux and the flux per pole are calculated, 
allowing subroutine MAGNET to compute the stator yoke and stator tooth flux densities, 
as well as ATTOT, the total magnetomotive force. This permits the calculation of a 
new, more accurate value of core loss WO and resistance RO as follows: 


WO = 3. * WFE* WSYOKE* ^*2 + WSTOTH* 

/BST\ * :|c «>n 

(30) 

\BK / 

Vbk/ J 

(3.*V2*V2) 


(31) 


WO 


where 


WSYOKE 

stator back-iron weight 

WSTOTH 

stator tooth weight 

BSY 

back- iron flux density 

BST 

stator tooth flux density 

V2 

airgap voltage (fig. 3) 


The process is then repeated using the latest value of RO until RO converges. 

When RO has converged, the ATTOT value calculated by MAGNET and the values 
of V2 and IMAG calculated by CIRCT are used to compute a new value of XO as fol- 
lows: 


XO 


V2 

0. 5* (IMAG + IMAG2) 


(32) 


where 


IMAG2 = (2. 22* P* ATTOT) 

(3.*N*KPS*KDS) 

With this new value of XO, the calculations for RO are then repeated. This double 
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convergence procedure is performed until both RO and XO have converged to their 
final values. These values are then held constant throughout the remaining motor an- 
alysis for all values of slip. 

This convergence procedure guarantees that the magnetizing ciirrent and the no-load 
core loss as obtained from the equivalent circuit analysis are the same as those calcu- 
lated from the magnetic material properties. 

In addition to computing the values of XO and RO, the double convergence proce- 
dure yields all no-load flux densities, no-load magnetomotive force across all parts of 
the magnetic circuit, no-load core loss, no-load magnetizing current, total useful flux, 
and total useful flux per pole. (Fig. 5 illustrates the no-load magnetic calculations in 
flow-chart format. ) 

Wind age loss calculations . - The windage loss calculations may be divided into two 
steps. The first step is to obtain the value of FWl, the windage loss at synchronous 
rotor speed. The second step, which is performed in subroutine CIRCT, is to calculate 
the windage loss at any other rotor speed from the equation 

FW = FWl* ( * C (33) 

VNSYNCH/ 


where 


FWl 

windage loss at S 3 mchronous speed 

FW 

windage loss at RPM 

RPM 

rotor speed 

NSYNCH 

synchronous rotor speed 

C 

constant (in the program, C =2.5) 


The value of FWl is allowable input to the program; if known, it may be read in 
directly. If FWl is not known, it may be omitted from the input data, in which case 
the program assumes it to be zero. Or, if a value of windage loss WL for a similar 
motor is known, WL may be read in. WL is called the reference windage loss; the 
motor for which the windage loss WL is known is called the reference motor. Assum- 
ing that sufficient additional information is provided to the program, the value WL will 
be scaled (multiplied by dimensionless parameters) to obtain the value FWl. The addi- 
tional information required by the program to scale WL is a complete description of 
the conditions under which the value WL was obtained. These conditions, called the 
reference conditions, consist of the following: 

(1) DIAREF, the diameter of the reference motor 

(2) LREF, the rotor length of the reference motor 
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(3) RPMREF, the rotor speed for which the value WL was obtained 

(4) GAPREF, the radial gap of the reference motor 

In addition, though not required, the reference viscosity VSCREF and the reference 
pressure PREF of the fluid in the reference motor cavity may be supplied. Then, 
assuming that the viscosity and pressure of the fluid in the motor being analyzed are 
also specified, further scaling with regard to these parameters will also take place. 

The exact format of the windage data required by the program is described in appendix A. 

Figure 6 is a flow chart that shows precisely how the value of FWl is determined 
by the program. Note that the program initializes the variables that specify the refer- 
ence conditions prior to reading the windage data. This initialization makes certain 
that all variables have numeric values, and it allows determination of which variables 
were omitted from the windage data. After having read the windage data the program 
checks if VSCREF was read in. If it was not, the program attempts to calculate a 
value from the constants CO to C4 and the temperature TREF. If CO to C4 were 
also omitted from the data, the value of VSCREF will remain zero and no scaling with 
regard to viscosity will be performed. 

The windage calculations described to this point are executed only once for each 
data set. By contrast, the remaining calculations are carried out once for each motor 
design data deck. 

Prior to reading the motor design data deck, a number of additional variables are 
initialized as shown in the flow chart. The design deck is then read and a heading 
"WINDAGE" is written on the output record. The program then checks if FWl was 
read in. If it was, its numerical value is printed out and all other windage calculations 
are bypassed. If FWl was not supplied, the program checks the value of WL. If WL 
equals 0, the program leaves FWl set to zero and proceeds as above. If WL is not 
zero, the program checks that numeric values for DIAREF, LREF, RPMREF, and 
GAPREF are supplied. If not, an error message is issued and FWl is left unchanged 
at zero. Otherwise, the program proceeds to compute 


FWl =WL* 


DR \ 
^DIAREF/ 


3.25 

* 


L Y /nSYNCh Y' / GAPREF ^- 
LREF/ \ NREF / \ G / 


(34) 


where 


DR 

NSYNCH 

L 

G 


rotor diameter of the motor being analyzed 
synchronous shaft speed, rpm 
rotor length 

radial gap between rotor and stator 
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Following this basic scaling operation the program will, if possible, scale FWl with 
respect to fluid viscosity and fluid pressure. This concludes the computation of FWl, 
the windage loss at synchronous rotor speed. The result, together with all pertinent 
data, is printed out. 

Calculation of windage loss at all rotor speeds other than synchronous is left to 
subroutine CIRCT using equation (32). 

Motor performance calculations . - Motor performance is computed at a given line- 
to-neutral voltage and line frequency. It consists of computii^ line current, shaft tor- 
que, output power, efficiency, power factor, input power, and internal losses. These 
performance parameters are computed by analyzing the motor’s equivalent circuit. 
Analyzing the circuit is equivalent to analyzing the motor. The performance parameters 
obtained directly from the equivalent circuit are line current, power factor, and input 
power. 

The power dissipated in the resistive element R1 is the ohmic loss of one phase of 
the armature winding. The power dissipated in RO is one-third of the motor core loss. 

The power W2 is three times the power dissipated in the resistive element 
(R2* 100)/S. It is given by 

W2 = 3. *12* 12* R2* 100. ^35^ 

S 

The symbol W2 represents both the loss in the rotor winding and the power delivered 
to the shaft. It is divided into its respective components as follows: 

W2 = 3. *12*12* R2 + 3. *12*12* R2* (100. - S) 

S 


where 

3. *12*I2*R2 power dissipated in the rotor windir^ 

3. *12*12* R2* (100. - S)/S gross shaft power 

Windage loss FW is subtracted from the gross shaft power to obtain the net power out- 
put from the motor at this point. The net shaft torque is then computed from the net 
power output and the shaft speed. 

All circuit analyses, as well as calculations of output power, efficiency, and losses 
are performed in subroutine CIRCT for one value of slip. The main program ESTDMTR 
increments the value of slip by 1 percent from zero to 100 percent. For each value of 
slip, subroutine CIRCT is called to perform all calculations described in this section. 
Following each call to subroutine CIRCT the main program prints the results of the cal- 
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culations for the one value of slip. 

If input data specify a value for rated torque, the naain program will check the com- 
puted value of torque after each call to subroutine CIRCT but before printing the results 
of the calculations. If the computed torque exceeds the specified rated torque, printout 
of the results is temporarily suppressed and normal processing is interrupted. The 
value of slip is no longer incremented by 1 percent. Instead the prc^ram searches, 
through an iteration process, for the value of slip that produces the rated output torque. 
Motor performance at that value of slip is then printed out. This output line is both 
preceded and followed by a blank line to offset it from the other output. Normal process- 
ing is then resumed with the value of slip at which processing was interrupted. Further 
comparison of computed torque with rated torque is discontinued. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 30, 1975, 

506-23. 
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APPENDIX A 


INPUT REQUIREMENTS 
Input Data Requirements 

The use of this computer program for the analysis of induction motors requires that 
the complete electromagnetic design be known. This includes physical dimensions, 
armature and rotor winding parameters, winding temperatures, and the magnetic char- 
acteristics of the materials to be used in the stator and the rotor. The design informa- 
tion is then transferred onto data cards for use with the program. A typical data deck 
is shown in figure 7. It consists of one or more data sets, all similar in makeup. Each 
data set contains, in the order required, two material data decks, windage data, and any 
number of motor design data decks. The material data decks must be in the order shown 
in the figure, that is, stator material data decks followed by the rotor material data 
deck. There must be two material data decks even if the rotor and stator are made of 
the same material. 

The material data decks and the windage data apply to all motor design data decks 
that follow, until a new material data deck is encountered. The appearance of a new 
material data deck within the input data signals the beginnii^ of a new data set. 


Preparation of Material Data Decks 

A material deck consists of from 5 to 16 cards as shown in figure 8. The first card 
is the material deck identification card. Its main purpose is to give the material name, 
which serves two functions: it identifies the material deck, and it is read by the com- 
puter and stored for later printout on the output record. The next four cards are the 
saturation curve data cards. They contain the coordinate values of as many as 14 arbi- 
trary data points located on the magnetization curve of the material. Following this are 
as many as 10 core-loss data cards - one for each lamination thickness at which core- 
loss calculations are anticipated. At least one core-loss data card is required for the 
stator material deck. The last card identifies the end of the material deck. Core-loss 
data cards and the last card are not needed for the rotor material deck. 

Each type of card that goes into the material data deck is described in detail here. 

Materi al deck identification card . - The material deck identification card must con- 
tain an *'M"-punch in card column 1. Column 2 should be blank. The material name 
should start in column 3 and may extend through column 80. 

Saturation curve data cards . - The fovm saturation curve data cards contain points 
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on the magnetization curve for the material. Each card is divided into eight 10-column 
fields. The entry in the first field of the first card must be the highest value of flux 
density of the selected points on the saturation curve. This is followed by paired values 
of magnetic flux density and magnetizing force. The values of flux density must be in 
ascending order. If less than 14 points are chosen for input such that less than fovir 
cards are required, blank cards must be inserted. The units must be kiloUnes per 
square inch for flux density and ampere-turns per inch for magnetizing force. During 
program execution, the original magnetization curve is approximately reconstructed by 
interpolation between points. The interpolation assumes a straight line on semilog paper 
between data points. 

Core- loss data cards . - The core-loss data cards are identified by the NAMELIST 
name FELOSS. Thus, they must contain the entry $FELOSS starting in card column 2. 
Other required entries on each card are of the form 

Variable name = Numeric value 

The permissible variable names and their definitions are shown in table I. The entries 
may be in any order, separated by commas. Each core- loss data card must be termi- 
nated by a "$ "-punch. 

The last core- loss data card must be 

$FELOSS LAST=.TRUE. $ 

Typical example . - Preparation of a material data deck is illustrated for M-19 
silicon steel. The first card of this material deck is shown in figure 9. 

To prepare the next four cards of the material data deck, the magnetization curve 
of the material is needed. The magnetization curve for M-19 steel is shown in figure 10. 
The 14 selected points are indicated by data symbols. The numeric values of these 
points are listed in the table insert of figure 10. The sequence in which the numbers are 
punched onto data cards is as follows: 116., 26., 1.30, 30., 1.45, ..., 110., 130., 
116., 135.0. 

Core loss as a function of frequency for M-19 steel is shown in figure 11. The fig- 
ure is for a lamination thickness of 0. 014 inch and a flux density of 64. 5 kilolines per 
square inch. The slope obtained from the figure is 

SLOPE = .^- (WVW2) 

In (F1/F2) 

where Wl, W2, FI, and F2 are arbitrarily selected points on the curve. Substituting 
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numeric values gives 


SLOPE = ^ (70. /I. 3) _ ^ 

In (1500. /lOO.) 

Thus, the core-loss data card for M-19 steel Avith a thickness of 0. 014 inch would ap- 
pear as 

$FELOSS WCORE=9.4, FCORE=400, BK=64.5, LT=0. 014, SLOPE=1.47 $ 

The value of SLOPE, plus one point on the core- loss-against- frequency curve defined 
by the coordinates WCORE and FCORE, allows the approximate reconstruction of the 
curve. Figure 11 shows how the program reconstructs the core-loss-against-frequency 
curve from these data. 

As many as nine additional cards for other lamination thicknesses may be added. 
The complete material deck for M-19 steel is shown in figure 9. 


Preparation of Windage Loss Data Cards 

A data card referencing the NAMELIST name WNDAGE must be included with each 
data set even if no windage calculations are to be performed by the program. If windage 
loss is to be neglected, no entry is required. If a known windage loss WL is to be scaled 
from the reference conditions to give the windage loss of the motor design at synchro- 
nous speed, the entries as defined in table II should be included. If it is desired to read 
in a value of windage loss at synchronous speed and b 3 q>ass the internal windage loss 
scaling operation, the entry FWl of NAMELIST name RATING (see next section and 
table III) should be used. For a description of windage loss calculations, see the main- 
text section Windage loss calculations. 


Preparation of Motor Design Data Deck 

The motor design data deck contains all the dimensions, the geometric configura- 
tion (in numerical code), and the winding parameters needed for an electromagnetic 
analysis of the motor design. The first card of the motor design deck is a title card 
similar to the first card of the material data deck. The first two columns of the title 
card must be blank. Any t 3 q>e of descriptive or identif 3 dng information may be punched 
in columns 3 to 80. This information is written on the output record for identification 
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purposes. The remaining cards of the motor design data deck are read with a READ 
statement referencing a NAMELIST name. For each NAMELIST name, one or more 
data cards are required to numerically define the variables included in that NAMELIST 
name. There are eight NAMELIST names; each is suggestive of the type of variable 
included in its list. Detailed information about each NAMELIST name is provided in 
table in, which lists all variables used with the alternator design data deck. All vari- 
ables belonging to the same NAMELIST name are grouped together. The NAMELIST 
names are arranged in the order in which the data cards must appear in the data deck. 
Units are given where applicable, and each variable is classified as either required or 
optional. Optional variables are read in at the discretion of the program user. In each 
case where an optional variable is omitted from the program input, an assumption re- 
garding that variable is made internal to the program. The remarks column of table III 
supplies specific explamtions in regard to all optional variables. 

To further clarify the definition of some variables, figures 12 and 13 are given. 
These figures are referenced in the table where applicable. 

A t 3 pical data card from the motor design deck is 

$STATOR D=1.07, L=1.24, SFS=0.909, DOS=2.50, LTS=0. 006 $ 

The card is for NAMELIST name STATOR and must be the second card following the 
motor design deck title card. 
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TABLE I. - ENTRIES ON CORE- LOSS DATA CARDS 


[namelist name, FELOSS.] 


Classification 

FORTRAN 

symbol 

Description 

Remarks 

Required^ 

WCORE 

Core loss at frequency FCORE, at flux 
density BK, and for lamination thick- 
ness LT, W/lb 



FCORE 

Frequency at which WCORE is deter- 
mined, Hz 

FCORE should be close to motor design 
frequency F to minimize the error of 
extrapolating WCORE from frequency 
FCORE to frequency F. 


BK 

Flux density at which WCORE is 
2 

determined, kilolines/in. 



LT 

Lamination thickness at which 
WCORE is determined, in. 


Optional 

LAST 

LAST=. TRUE . indicates end of 
core -loss data set. 

This variable is required on the last 
card of core- loss data set. It must be 
the only variable on that data card and 
it must be set equal to . TRUE. If this 
variable appears on any other card, it 
must be set equal to . FALSE. 

! 

SLOPE 

Slope of core loss-against-frequency 
plot (assumed linear on log- log paper) 
for constant flux density BK; if W1 
is core loss at frequency FI and W2 
is core loss at frequency F2, 

Wi 
In— t 

Wo 

SLOPE = ^ 

In— i 

F2 

If SLOPE is omitted from input, the 
program will assume a value of slope 
given by 

SLOPE = 1- t.l64.*LT 
1. + 82.*LT 

SLOPE may also be omitted if F of 
NAMELIST name RATING equals 
FCORE. 



(fig. 11) 



^For stator material only. 
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II 


TABLE n. - ENTRIES ON WINDAGE-LOSS DATA CARDS 
[namelist name; WNDAGE.] 

FORTRAN Description Remarks 

symbol 

WL Windage loss at reference conditions, If the windage loss for the motor being 

W analyzed is known, its value may be 

DIAREF Reference diameter, in. read in by using the variable FWl of 

LREF Reference length, in. NAMELIST name RATING (table III). 

RPMREF Reference rotational speed, rpm If FWl is used, the scaling will be by- 

GAPREF Reference gap length, in. passed and hence no entries on the 

WNDAGE data card are required. 

VSCREF Reference viscosity, Ibm/ft-sec If VSCREF is omitted, the program will 

attempt to calculate a value based on the 
polynomial 

VSCREF = CO + C1*TREF 

+ C2*TREF**2 + C3*TREF**3 
+ C4*TREF**4 

If the result of the calculation is <0, no 
windage loss scaling with respect to vis- 
cosity will be made. 

CO These constants are coefficients of a Any or all constants omitted from the 

Cl fourth- degree polynomial that gives input will be assumed to equal zero. 

C2 the viscosity as a function of temper- 

C3 ature. The constants must be chosen 

C4 such that viscosity is in Ibm/ft-sec for 

temperature in °C. 

TREF Fluid temperature at reference condi- TREF is the temperature for which 

tions, °C VSCREF is calculated if necessary. 

PREF Pressure of fluid in airgap at reference If omitted, scaling of windage loss with 
conditions, Ib/in. pressure is not possible. 

^These variables are required only if it is desired to scale the windage loss from a known condition to the 
motor being analyzed. If scaling is not desired, no entries are required. 


Classification 

Required^ 


Optional 
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TABLE m. - INPUT REQUIREMENTS FOR MOTOR DESIGN DECK 


NAMELIST 

Classification 

FORTRAN 

Description 

Remarks 

name 


symbol 



RATING^ 

Required 

NSYNCH 

Synchronous speed, rpm 




F 

Line frequency, Hz 




VI 

Line -to- neutral voltage, V 



Optional 

TRATED 

Rated torque, in. - lb 

Normally motor characteristics are 
computed at predetermined values of 
slip only. If a value of TRATED is 
specified, motor characteristics at the 
value of slip corresponding to torque 
TRATED will also be computed. 



FWl 

Windage loss at synchronous 

If specified, the internal windage loss 




speed, W 

scaling will be bypassed. 



xo 

Reactance values of induction 

If any of these reactances and resist- 



XI 

motor equivalent circuit 

ances are specified, internal calcula- 



X2 

(fig. 3), G 

tions for that circuit element are by- 



RO 

Resistance values of induction 

passed. If all are specified, all 



R1 

motor equivalent circuit 

internal calculations are b 5 q)assed ex- 



R2 

(fig. 3), G 

cept for the equivalent circuit analysis. 
In this case, all remaining data cards 
should be omitted. (The material data 
decks and the $WNDAGE data set are 
still required; the data therein are not 





used. ) 

STATOR 

Required 

D 

Stator inside diameter, in. 




L 

Stator stack length, in. 




LTS 

Stator lamination thickness, 



1 

i 

I 

DOS 

in. 

Stator outside diameter, in. 



1 

Optional 

SFS 

Stator stacking factor 

If omitted, stacking factor is calculated 
as follows: 





SFC — LTS 

~ LTS + 0. 0005 

SSLOTS 

Required 

SSTYPE 

Stator slot type (Choose from 
t 5 rpe 1 to t5T>e 6, as shown in 
fig. 12.) 

SSTYPE is required for all stator slots. 


^This card must be preceded by a title card (see input, p. 43. ) 
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NAMELIST Classification FORTRAN 
name symbol 


TABLE ni. - Continued. 
Description 


SSLOTS 


symbol 

Required DIS 


DIS 

or 

DSS 

or 

SCAREA 

or 

CSRATO 

Conductor depth (fig. 12), in. 
Slot depth (fig. 12), in. 

Slot area (fig. 12), in. ^ 

Space factor 

D6S 

Stator slot dimension (fig. 12), 
in. 

QS 

Number of stator slots 

D2S 

DSS 

WSS 

WSS6 

Slot dimension (fig. 12(a)) 

D2S 

DSS 

WSS6 

WSS 

or 

Slot dimension (fig. 12(b)) 

STWDTH 

Stator tooth width if constant, 


Remarks 


One of these variables is required for 
all stator slot types. 


D6S and QS are required for all stator 
slots. 


One of these variables is required for 
stator slot type 1 only. 


Required for stator slot type 2 only. 


Slot dimension (fig. 12(c)), in. Required for stator slot type 3 only. 


D2S 

D3S 

D4S 

D5S 

WSSl 

WSS6 

WSS2 

or 

STWDTH 


Slot dimension (fig. 12(d)), in. Required for stator slot type 4 only. 


Stator tooth width if constant. 


TABLE ni. - Continued. 


NAMELIST 

name 

Classification 

FORTRAN 

symbol 

Description 

Remarks 

SSLOTS 

Required 

D4S 

WSSl 

Slot dimension (fig. 12(e)), in. 
Slot dimension (fig. 12(e)), in. 

Required for stator slot tjqae 5 only. 



D2S 

D3S 

D4S 

D5S 

WSSl 

WSS6 

WSS2 

or 

STWDTH 

Stator slot dimension 
(fig. 12(f)), in. 

Stator tooth width, in. 

Required for stator slot type 6 only. 


Optional 

PfflS 

One- half of angle at which 
sides of stator slots diverge 
(fig. 12), deg 

PHIS is allowable input only if WSS2 is 
specified and STWDTH is omitted for 
slot types 2, 4, or 6. If, in that case, 
PHIS is also omitted, constant tooth 
width is assumed. 



WSS2 

Stator slot dimension 
(fig. 12(e)), in. 

Optional for slot type 5 only 

STRWDG 

Required 

CSS 

PC 

B 

SPITCH 

ASTRND 

or 

AWG 

S 

Number of conductors per 
stator slot 

Number of parallel circuits 
Armature coil extension 
(fig. 12), in. 

Stator winding pitch, per unit 

Cross-sectional area of stator 
2 

conductor strand, in. 

Gage size of stator conductor 
strand, AWG 

Clearance between armature 
coils at end turns (fig. 13), in. 



Optional 

SWMAT 

Stator winding material code 

SWMAT = 1 for aluminum 

SWMAT = 2 for brass 

SWMAT = 3 for copper 

Unless one of the other two is specified, 

SWMAT = 3 is assumed. 



LS 

Length of one armature con- 
ductor (one- half of armature 
coil length), in. 

If omitted from input data, program 
calculates value internally. The cal- 
culations assume a form -wound 
armature. 
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TABLE in. - Continued. 


NAMELIST 

name 

Classification 

FORTRAN 

symbol 

Description 

Remarks 

STRWDG 

Optional 

TSW 

Temperature of stator wind- 
ing, °C 

A temperature of 25° C is assumed un- 
less otherwise specified. 

STRNDS 

Number of strands per arma- 
ture conductor 

One strand per conductor is assumed. 

ROTOR 

Required 

LTR 

DIR 

Rotor lamination thickness, 
in. 

Rotor lamination inside 
diameter, in. 

If omitted from input data, program 
assumes it to equal rotor slot pitch or 
stator slot pitch, whichever is greater. 

If omitted, stacking factor is calcu- 
lated as follows: 

SFR ~ LTR 

LTR + 0. 0005 

Optional 

SKEW 

Rotor slot skew, in. 

SFR 

Rotor lamination stacking 
factor 

RSLOTS 

Required 

RSTYPE 

Rotor slot type (choose from 
type 1 to type 6 as shown in 
fig. 12) 

Required for all rotor slots. 

SB 

or 

DIR 

or y 

DSR j 

Rotor bar cross-sectional area 

(fig. 12), in. 2 

Slot depth (fig, 12), in. 

One of these is required for all rotor 
slots 

Required for all rotor slots 


D6R 

Slot dimension (fig. 12), in. 

D2R 

WSR 

WSR6 

Slot dimension (fig. 12(a)) 
Slot dimension (fig. 12(a)) 
Slot dimension (fig. 12(a)) 

Required for rotor slot type 1 only. 

D2R 

WSR6 

WSR 

or 

RTWDTH 

Slot dimension (fig. 12(b)) 
Slot dimension (fig. 12(b)) 
Slot dimension (fig. 12(b)) 

Rotor tooth width if constant, 
in. 

Required for rotor slot type 2 only. 
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TABLE m. - Continued. 


NAMELIST 

name 

RSLOTS 


RTRWDG 


Classification 


Required 


Optional 


Required 


FORTRAN 

Description 

Remarks 

symbol 




D2R 

Slot dimension (fig. 12(c)), in. 

Required for rotor slot type 3 only. 

D3R 

D4R 

WSR 

WSRl 

WSR6 




D2R 

Slot dimension (fig. 12(d)), in. 

Required for rotor slot type 4 only. 

D3R 




D4R 




WSRl 




WSR6 




WSR2 




or 




RTWDTH 

Rotor tooth width if constant, 



in. 



D4R 

Slot dimension (fig. 12(e)), in. 

Required for rotor slot type 5 only. 

WSRl 

Slot dimension (fig. 12(e)), in. 


D2R 

Rotor slot dimension 

Required for rotor slot type 6 only. 

D3R 

(fig. 12(f)), in. 


D4R 




WSRl 




WSR6 




WSR2 




or 




RTWDTH 

Rotor tooth width, in. 


WSR2 

Rotor slot dimension 

H omitted, WSR2 = DSR - D4R for 


(fig. 12(e)), in. 

slot type 5 only. 

PHIR 

One -half of angle at which 

Allowable input only if RTWDTH is 


sides of rotor slots diverge 

omitted for slot type 2, 4, or 6 and 


(fig. 12), deg 


WSR2 is specified. If, in that case, 
PHIR is also omitted, constant tooth 
width is assumed. 

NB 

Number of rotor bars 


TER 

End- ring thickness (measured 
in axial direction), in. 


BR 

Axial clearance between end- 
ring and rotor laminations, in. 
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TABLE m. - Concluded. 


NAMELIST 

name 

Classification 

FORTRAN 

symbol 

Description 

Remarks 

RTRWDG 

Optional 

LB 

Rotor bar length, in. 

If omitted, LB is calculated as follows: 





LB = + SKEW^ + 2* (BR + TR) 



DERI 

End-ring outside diameter, 
in. 

If omitted, DERI is calculated as fol- 
lows: 





DERI = DR - 2. * (D4R + D3R) 



DER2 

End-ring inside diameter, in. 

If omitted, DER2 is calculated as fol- 
lows: DER2 = DR - 3.*DSR or 
DER2 = 1. 1*DIR whichever is greater. 



RWMAT 

Rotor winding material code 

RWMAT = 1 for aluminum 
= 2 for brass 
= 3 for copper 

Unless one of the other two is specified, 
RWMAT = 3 is assumed. 



TRW 

Rotor winding temperature, ®C 

A temperature of 25° C is assumed un- 
less otherwise specified. 

AIRGAP 

Required 

G 

Airgap, in. 



Optional 

TFLUID 

Temperature of fluid in air- 
gap, °C 

If omitted, it is assumed that 
TFLUID = TREF. (See NAMELIST 
name WNDAGE.) 



VSCFLD 

Viscosity of fluid in motor 
cavity, Ibm/ft-sec 

If omitted, program will calculate value 
of VSCFLD based on temperature 
TFLUID and the constants CO to C4 
of NAMELIST name WNDAGE. If all 
values CO to C4 are omitted from 
the input data, the results of the calcu- 
lation will be VSCFLD = 0. In this 
case, the windage loss WL will not be 
scaled with regard to viscosity. 



PFLUID 

Pressure of fluid in airgap, 
psi 

This needs to be included only if it is 
desired to scale the windage loss value 
WL to the new pressure level PFLUID. 



FLDNME 

Name of fluid in airgap (may 
be a maximum of six charac- 
ters long and must be enclosed 
in single quotation marks) 

If specified, the name will be printed 
on the output record. No other action 
occurs. 


28 












APPENDIX B 


TYPICAL INPUT AND RESULTANT OUTPUT 

A complete data set identifying the material deck, the windage data, and the motor 
design deck is given here: 


PvRMfNDUR 


1S4 . 

1?.9 

1 .92 

38.7 2 

.62 

77. 

4 3.23 

90.3 


3-53 

103. 

4 . 38 

109,7 S 

.25 

116. 

^ *66 

1 22 . 

5 

S. 6 P 

129 . 

12 .8 

1 55.5 

20.2 

142. 

44.4 

145. 

3 

10 1 . 

146.5 

363. 

184 . 

202C. 





SFELCSS 

i<cFPE-?i -n , 

FC 0 RF= 8 ;: 2 . 

, SL CPF -1.22 

, SK=77.4 

, LT 

= 0.006 5 



5FLL 0*^S 

WCcPE-34.5, 

KCORFzFO 3. , 

» SL0PE=1.34 

, BK=77.4 

• LT 

=o.no8 s 



SF EL Ocs 

wCORF- 3 T .rj , 

F C 0 Rf r s D n . , 

r SL0 Pf=1.4S 

, PKz77.4 

, lt 

rO.OlO S 



SF LLOSS 

UG ore =42,0, 

F C 0 R E' = 9 iJ , 

. SL0PE=1.57 

, BK-77.4 

. LT 

=0.014 ^ 



SFE L C'^S 

L ^ STr . TkUt 

. -i. 







PSkHENL-Uk 








I*- 4 . 

12.9 

1 .92 

GH.T 2 

.62 

77. 

4 3.23 

90.3 


3. F 3 

i:)3 . 

4 . 3S 

129.7 5 

.25 

11 b. 

6.66 

122 . 

C 

F . 6 P 

129. 

1 2 .F 

1 3S .5 

?D- 2 

142. 

44.4 

145. 

3 


U'l. 148.3 1S4. r320, 

4FELCSS kvCf-RFzZl.n, FCORF , SLOPE -1.22, Bk-77.4, LT= 0 .006 1 

SFELO'^S hLC^^'E-2‘^ » FCOFf = tO':., SloPF = 1-34, BK-77.4, LTr[T.Q03 % 

SFCLCSS WCoOL-30.3, FCORH-rjO., SL0PLri.q5, BK-77.4, LTm.OlO % 

5FLL0SS 'r^C CRE -4C .0 , F L f RE - 8 [) 3 . , SLOPE-l.S?, BK-77.4, LT-D.CIM 1 

IFELOSS LiSTr.TPUF. 

SWhOflrE WL=45.0, DIflRlF=l.:s, LRFE=1.126, PPHRlF=120Ql, GAPREFrO.ClO, 
Cn=1.728E-3, Cir-3 .r2^>3:’l -3 , C 2 1: 0 . 3 8 7 6 E -6 , C 3--0 . 3294E -fi , C 4 -0 . 1 08 E -1 C , 
TREF z 2o. 

12fHJ H2 COOLANT FU^'P ‘'OTOR 

4WATIRG N SY^CH^ 1 2200 , F::l?ir, ,/ 1 - 1 2 ' , TRATEH^Z.u % 

SST^TOR :=1.C7, L-\.Z^t SrS=3.909» OOSzZ-^u, LTS=2.D06 * 

SSSLUTS C^"30, SSTYP£ = 6» 02S:: ' .1.:, D4S = L.. D.iSrC., 

KSSlrr:,04^, 06S-G-G1C8, S T u F T Kr 1 . r q G , OSS-G-O^G, 

U6S = j.:iS, USS6Z''’ . J 1 '■< 3 

SSTUU'DG A^C = 25, CSSr56, PC-^, ^ . 1 , S'-KAT^i., S P I T C H - . 66 6 7 , S-O.OIOS, 

TSWrir^., LS = 2.03 % 

$RGTOP lTr-^.Q26, S'"P-?.9i, Dir “ , 14 F ^ 

SR SLOTS I- S TYPE = 1, D4r?zJ.J50, FiR^r., WSF'lrQ.LHS, r>2R = 0 .C 0 S, 

LSRra.^<)2, WSR:^.nS3. 06 Rz.j.» WSk6-''.(?J3 S 

SRTPWno Dc-29, TER-.13S, RWMAT r?., FRr^., OrRi=l.ni3, DERprQ.SO, TRW-3D. 

SAJRGAp r- = ^.GC6, TFLUlO r 20., FL^J^»Mfr •DC-200* % 



Stator 

material 

data 


Rotor 

material 

data 


Windage data 
Title card 


Motor 

design 

data 
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The output that resulted from using this data set with the induction motor computer pro- 
gram is as follows: 


STATOR MATERIAL 

VANADllH PERHENDUR 


IKILOLINES/SO-INI 


f A-TURN/INl 


12.90 

38.70 

77.i|0 

90-30 

103.00 
109.70 

116.00 
122 .50 
129.00 
135.50 
1 <| 2.00 
IMS. 30 
1M8.30 
15m .00 


1.92 
2.62 
3.23 
3-53 
M.35 
5.25 
6.66 
6.68 
12-50 
20.20 
MM. MO 
101.00 
363.0a 
2020.00 


CORE-LCSS DATA 
CORE-LOSS 

LAM THK 

FREO 

FLUX DNSTY 

SLOPE 

21.0 

.006 

800.0 

77.4 

1.2 

2M.5 

• 0 Q 8 

800 .0 

77-4 

1.3 

30.0 

.010 

800.0 

77.4 

*1.5 

40.0 

.014 

800 .0 

77-4 

1.6 


ROTOR MATERIAL 

VANADltM PERMENDUR 


IKILOLINES/SO-INI 


(A-TURN/INI 


12.90 

1.92 

38 -70 

2.62 

77.40 

3-23 

90 .30 

3.53 

1C3.00 

4.35 

1C9 .70 

5.25 

116.00 

6.66 

122 .50 

8.68 

129 .00 

12.50 

135 .50 

20.20 

142 -00 

44.40 

145-30 

101.00 

148 .30 

363.00 

154-00 

2020.00 


1200 H2 COOLANT PUHP MOTOR 


RATING 


SYNCHRONOUS SPEED 

12000. RPM 

FRECUENCY 

1200. HZ 

POLES 

12. 

L-N VOLTAGE 

120.0 VOL 

RATED TORQUE 

2.0 IN-I 

TOR 

BORE diameter 

1 .070 

OtTSlDE diameter 

2.500 

DEPTH BELOy SLOT 

.195 

LENGTH 

1 .240 

lamination THICKNESS 

.006 

STACKING FACTOR 

.909 

STATOR IRON WEIGHT 

.743 
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STA10R SLOTS 


SLOT TYPE 6 

NO. OF slots 

36. 

TOOTM WIDTH .045 

SLOT DEPTH 

.520 

ySSl .048 

DIS 

.405 

USS2 .048 

D2S 

.100 

USS3 .128 

D3S 

.ODD 

USS4 .066 

D4S 

*000 

USS5 .128 

OSS 

.010 

yss6 .010 

D6S 

.D15 

USABLE AREA .029 

TOTAL AREA 

.047 

SPACE FACTOR .491 



STATCR UINDING 



MATERIAL 


3 

CONDUCTORS PER SLOT 


56. 

PARALLEL CIRCUITS 


2. 

PITCH 


.667 

AXIAL EXTENSION BEYOND 

CORE 

.100 

conductor CROSS-SECTIOn 


.252 

STRAND CROSS-SECTION 


.252 

CONDUCTOR LENGTH 


2.530 

CLRNCE BTWN END-TURNS 


.010 

temperature «c> 


30. 

AXIAL END-TURN LENGTH 


.730 

OVERALL ARMATURE LENGTH 


2.700 

PITCH FACTOR 


.866 

DISTRIBUTION FACTOR 


1 .000 

armature WEIGHT 


.412 

TOTAL armature WIRE LENGTH 

425.040 

strands/conductor 


1 . 

strand SIZE 


25 

ROTOR 



ROTOR diameter 


1.058 

INSIDE DIAMETER 


.450 

LAMINATION THICKNESS 


.006 

STACKING Factor 


.910 

SLOT SKEW 


.115 

CEPTH BELOW SLOT 


.212 

ROTOR IRON weight 


.185 

ROTOR SLOTS 



SLOT type 1 

NO* OF SLOTS 

29. 

SLOT WIDTH .053 

SLOT DEPTH 

.092 

WSR] .000 

DIR 

.087 

WSR2 -000 

02R 

.005 

WSR3 .000 

D3r 

•DOO 

WsR4 .053 

D4R 

.000 

WSR5 .053 

D6R 

.000 

WSR6 .003 



USABLE AREA .004 

TOTAL area 

.005 

ROTOR winding 



material 


3 

BAR length 


1.515 

BAR CROSS-SECTION 


.004 

end-ring OUTSIDE DIA 


1 -013 

END-RING INSIDE DIA 


.500 

end-ring THICKNESS 


.135 

STACK-TO-ENO-RING CLRNCE 

.000 

WINDING temperature CC) 


30. 

WEIGHT 


.100 

COMPONENT OF R2 CUE TO 

BARS 

2.084 

COMPONENT OF R2 DUE TO 

END RINGS 

.071 





A1R64P 

ACTUAL AIRGAP 
EFFECTIVE AIRGAP 
HAGRETIZING REACTANCE 

LEAKAGE REACTANCES fOHHI 

SLOT 

END-CONNECTION 

SKEU 

ZIG-ZAG 

peripheral 

WEIGHT 

total (electrohagnetji 


•0060 

•omi 

CAIR GAP ONLYI 12.73 

STATOR ROTOR 

8.332 1.673 

.579 .161 

.877 .877 

.619 1.173 

.030 

) 1.440 


STATOR material - VANADIUM PERMENDuR 
B MAX □ 154. 

CORE LOSS at 77.4 KL/S0-1N= 34.4 W/LB 


ROTOR material -- VANADIUM PERMENDUR 
B MAX = 154. 

magnetization characteristics 
ino-load, rated voltage I 

total USEFUL FLUX 158.42 KILOLINES 

USEFUL FLUX/POLE 8.41 


FLUX DENSITIES 
AIRGAP 

STATOR TOOTH 
STATOR YOKE 

rotor tooth 

ROTOR YOKE 


38.01 kL/SO-IN 

86.76 

19.13 

88.07 

17.58 


AMPERE-TURNS PER POLE 
AIRGAP 

STATOR TOOTH 
STATOR YOKE 
ROTOR TOOTH 
ROTOR YOKE 


167.53 

1.79 

•62 

.32 

.18 


TOTAL 


170.44 


MAGNETIZING CURRENT 
AIRGAP VOLTAGE 
N.L. CURRENT DENSITY 
CORE LOSS 


5.22 AMPERES 
65.19 
10380.04 

41. WATT 


WINDAGE 


DESIGN REFERENCE 

CONDITION CONDITION 


WINDAGE LOSS* W 56. 

DIAMETER 1.058 

Length i.24D 

RPM 12000. 

GAP .006 

TeHP, deg C 25. 

VSCSTY# LBM/FT-SEC .121-02 

PRESSURE, LBS/SO-IN .000 


45 . 

1.050 

1.125 

12Q00. 

.010 

20 . 

.129-02 
• 000 


FLUID 


DC-200 
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EOUIVALENT 

CIRCUIT parameters 



R1 = 

1.190 

XI r 

10.43b 

R2 : 

2.155 

X2 = 

3.08^ 

RO = 

312.278 

XD □ 

12.51u 


HOTOR performance at 120.00 VOLT, 1200. D HZ 


TOROUE 

SLIP 

RPH 

P 

-OUT 

1 

EFF 

PF 

P-IN 

PRi 

SEC 

IRON 

FW 

UN-LBS> 

{PERCENT ) 

IHP ) 

(WATT ) 

< amp } 

{Percent 1 1* 

=LEAOInGI (UATT) 

LOSS 

tWATT) 

LOSS 

(WATT » 

LOSS 

IWATIl 

(MATTI 

• 03 

1 .00 

11880.00 

.00 

3.69 

5 .23 

1.87 

.10 

197.23 

97.67 

.59 

40 .66 

54.62 

.44 

2.00 

11760.00 

.08 

61 .05 

5 .26 

23. P7 

• 14 

255.72 

98.80 

2.33 

40.29 

53.26 

• 84 

3.00 

11640. DO 

.15 

115.97 

5.31 

36.90 

.16 

312.96 

100.61 

5.18 

39.80 

51.91 

1.23 

4.00 

11520.00 

.22 

166.48 

5.37 

45.19 

.19 

368.36 

103.05 

9.04 

39.20 

50.58 

1 .55 

5.00 

114GC.00 

.29 

213.69 

5.45 

50.71 

.21 

421.39 

106.08 

13.84 

38.51 

49.27 

1.94 

6.00 

11280.00 

.34 

256.04 

5.5 4 

54 .46 

.24 

471 .63 

109.61 

19.46 

37.73 

47.99 

2.00 

6.20 

11255.47 

• 36 

265.15 

5 .56 

55 .07 

.24 

481.52 

110.39 

20.70 

37.56 

47.73 

2.25 

7.00 

11160.00 

.40 

295.77 

5 .64 

57.02 

.26 

518.73 

113.58 

25.78 

36.88 

46.72 

2.54 

8.00 

11040.00 

.44 

330*41 

5.75 

58.74 

.27 

562.46 

117.92 

32.69 

35.97 

45.47 

2.81 

9.00 

10920.00 

• 48 

360.81 

5.86 

59.87 

• 29 

602*57 

122.54 

40.06 

35.02 

44.25 

3.05 

10.00 

108C0.00 

.52 

387.05 

5.97 

6Q.54 

. 30 

639.30 

127.38 

47.79 

34 .04 

43.04 

3.26 

11.00 

1C680.Q0 

.55 

4Q9.32 

6.09 

60,88 

.31 

672.36 

132.37 

55.76 

33.05 

41.86 

3.44 

12.00 

10560.00 

• 57 

427.85 

6.21 

60 . 05 

.31 

701.93 

137.46 

63.89 

32.04 

40.69 

3.6l 

13.00 

10440.00 

.59 

442.86 

6.32 

60.82 

.32 

728.12 

142.58 

72.08 

31.04 

39.55 

3.74 

14.00 

10320.00 

.61 

454 .66 

6.43 

60 -^3 

.32 

751.09 

147.70 

80.27 

30.05 

38.42 

3.86 

15.00 

10200.00 

.62 

463.51 

6.54 

60.11 

.33 

771.04 

152.76 

86.38 

29 .08 

37.31 

3.96 

16.00 

1008D.D0 

.63 

469.70 

6.65 

59.59 

.33 

788.17 

157.75 

96.37 

28.13 

36.22 

4.04 

17.00 

9960.00 

.63 

473.52 

6 .75 

58 .99 

.33 

802 . 70 

162.63 

104.19 

27.21 

35.16 

4.10 

18.00 

9840.00 

.64 

475.22 

6 .85 

58.32 

.33 

814.83 

167.39 

111.80 

26.31 

34*11 

4.15 

19.00 

9720.00 

.64 

475. U7 

6.94 

57.80 

• 33 

824.79 

172.00 

119.19 

25 .45 

33.08 

4.19 

20.00 

9600.00 

.63 

473.^9 

7. 03 

56.83 

.33 

832.77 

176.46 

125.34 

24.62 

32.06 

4.21 

21.00 

94 80.00 

.63 

470.11 

7.12 

56.03 

.33 

838.99 

180.75 

133.22 

23.83 

31.07 

4.23 

22.00 

9360.00 

• 62 

465.72 

7.20 

55.21 

.33 

843.61 

184 .88 

139.85 

23.07 

3o*10 

4.23 

23.00 

9240.00 

.62 

460.30 

7.27 

54.36 

.32 

846.82 

188.84 

146.20 

22.34 

29.14 

4.23 

24.00 

9120.00 

• 61 

454 .01 

7.35 

53.49 

.32 

848.77 

192.63 

152.28 

21 .65 

28.21 

4.22 

25.00 

9000.00 

• 60 

446.99 

7.41 

52.61 

.32 

849.62 

196.26 

158.09 

20.99 

27.29 

4.21 

26.00 

8880.00 

.59 

439.3? 

7.48 

51.72 

.32 

849.49 

199.72 

163.65 

20.37 

25*39 

4.18 

27.00 

8760.00 

.58 

431.27 

7.54 

50.83 

.31 

848.51 

203.02 

168.94 

19.77 

25.50 

*♦.16 

28.00 

8640.00 

.57 

422.77 

7.60 

49.93 

.31 

646.78 

206*16 

173.99 

19.21 

24.64 

9.13 

29.00 

8520.00 

.55 

413.97 

7.65 

99.03 

.31 

644 .40 

209.16 

178.80 

18.67 

23.79 

4.10 

30.00 

8400.00 

.54 

4o4 .94 

7.71 

48.32 

.30 

841.47 

212.01 

183.39 

18.17 

22.96 

4.06 

31.00 

8280.00 

.53 

395.74 

7.76 

47,22 

.30 

836.05 

214.73 

187.75 

17.69 

22.15 

4 .Q2 

32.00 

8 1 60.00 

.52 

386.42 

7.80 

46.32 

.30 

834.22 

217.31 

191.90 

17.23 

21.36 

3.99 

33.00 

8040.00 

.51 

377.05 

7.85 

45.43 

.29 

830.05 

219.77 

195.85 

16.80 

20.58 

3.95 

34.00 

7920.00 

.49 

367.65 

7.89 

44 .53 

.29 

825.58 

222.11 

199.61 

16.39 

19.82 

3.90 

35.00 

7800.00 

• 48 

358.^6 

7.93 

43 .64 

.29 

820.86 

224.34 

203.18 

16.00 

19.08 

3.86 

36.00 

7680.00 

.47 

348 .92 

7.96 

42.76 

.28 

8l5.9q 

226*46 

206.59 

15.63 

18.35 

3.82 

37.00 

7560.00 

.46 

339 .63 

8 .00 

4 1.89 

• 28 

810.86 

228.47 

209.83 

15.28 

17.65 

3.77 

38.00 

7440.00 

.44 

330.44 

6.03 

41.02 

.28 

805.65 

230.39 

212.92 

14.94 

16.95 

3.73 

39.00 

7320.00 

.43 

321.34 

8.07 

40.15 

.28 

800-33 

232.22 

215.86 

14.63 

16.28 

3.69 

40.00 

7200.00 

.42 

312.37 

8.10 

39.29 

• 27 

794.94 

233.96 

218.66 

14.33 

15.62 

3.64 

41.00 

7080.00 

• 41 

3d3.52 

8.12 

38.44 

.27 

789.49 

235.62 

221.33 

14.04 

14.98 

3.6C 

42. DO 

6960.00 

.40 

294.81 

8.15 

37.60 

.27 

784.01 

237.21 

223.87 

13.77 

14.35 

3.56 

43.00 

6840.00 

.38 

286 .i:4 

8.18 

36.77 

.26 

778.51 

238.71 

226.30 

13.51 

13.74 

3.51 

44.00 

6720.00 

.37 

277.82 

8 .20 

35.94 

.26 

773.01 

240.15 

228.62 

13.27 

13.15 

3.47 

45.00 

6600 .00 

.36 

269.56 

8.23 

35.12 

.26 

767.52 

241.53 

230.83 

13.04 

12.57 



3 . M3 

96.00 

6980.00 

.35 

261 .45 

8.25 

39 .31 

• 26 

762.05 

292 .84 

232.94 

12.81 

12.00 

3.39 

97.00 

6360.00 

.39 

253.50 

8 . 2 ? 

33.50 

.25 

756.61 

294 .10 

234.96 

12.60 

11.95 

3.35 

98.00 

6290.00 

• 33 

295.71 

8.29 

32.71 

.25 

751.21 

295.29 

236.89 

12.90 

10.92 

3.31 

H9 .00 

6120.00 

.32 

238 . L7 

8.31 

31.92 

.25 

795 .86 

296 .99 

238.73 

12.21 

10.40 

3.27 

50 ,00 

6000.00 

.31 

230.60 

8.33 

31.19 

.25 

790.56 

297.59 

240.50 

12.03 

9 . 9Q 

3.23 

51.00 

5880.00 

.30 

223.28 

8.39 

30.36 

.29 

735.32 

298.59 

242.19 

11.85 

9.41 

3.19 

52.00 

5760.00 

.29 

216.11 

8.36 

29.60 

.29 

730.13 

249 .59 

243.81 

11 .68 

8.94 

3.15 

53.00 

5690.00 

.28 

209.10 

8.38 

28 .89 

.29 

725.01 

250.55 

245.36 

11.52 

8.48 

3.11 

54.00 

5520. DO 

.27 

202.23 

8.39 

2 8 . C9 

.29 

719.96 

251.48 

246.89 

11.37 

8.04 

3.08 

55.00 

5900.00 

.26 

195.52 

8.9 1 

27.35 

.29 

719.98 

252.36 

248.27 

11-22 

7.61 

3 . C4 

56.00 

5280.00 

.25 

188.95 

8.92 

26.61 

.23 

710.07 

253.21 

249.63 

11 .08 

7.19 

3.01 

57.00 

5160.00 

.29 

182.52 

8 .99 

25.88 

.23 

705.23 

254.02 

250.95 

10.95 

6.79 

2.97 

58.00 

5090.00 

.29 

176.23 

8.95 

25 .16 

.23 

700.97 

254.81 

252.21 

10.82 

6.40 

2 . 9R 

59.00 

9920.00 

.23 

170 . U8 

8 . 4b 

29.44 

.23 

695.78 

255.56 

253.42 

lo-'^o 

6-03 

2.90 

60.00 

4800.00 

.22 

169 . U5 

8.97 

23.79 

.23 

691.17 

256.28 

254.58 

10.58 

5.67 

2.87 

61.00 

9680.00 

.21 

158.16 

8.98 

23 . P3 

.?2 

68 6 *63 

256.97 

255.70 

10.97 

5.32 

2.89 

62.00 

9560.00 

.20 

152.39 

8.50 

22.39 

.22 

682.16 

257 .64 

256.78 

10.36 

4.99 

2. 81 

63.00 

9490. OD 

.20 

196.75 

8.51 

21.65 

.22 

677.77 

258.28 

257.81 

10-26 

4.66 

2.78 

69 .QO 

9320.00 

.19 

191 .23 

8.52 

20.97 

.22 

673.45 

258.90 

258 . 8 1 

10.16 

4.36 

2.75 

65.00 

9200.00 

.18 

135.82 

8.53 

20.30 

.22 

669.21 

259.50 

259.77 

10.06 

4.06 

2.72 

66.00 

9080.00 

.17 

130.52 

8 .54 

19.63 

.22 

665.04 

260 .07 

260.70 

9.97 

3.78 

2.69 

67.00 

3960.00 

.17 

125.39 

8 *54 

18.96 

.21 

660.94 

260.63 

261.59 

9-88 

3,50 

2.66 

68.00 

3890.00 

.16 

120.26 

8 .55 

18.31 

.21 

656.92 

261-16 

262.45 

9.80 

3.24 

2.63 

69.00 

3720.00 

.15 

115.29 

8.56 

17.66 

.21 

652.96 

261.68 

263.28 

9.71 

3 . Q0 

2.61 

70.00 

3600.00 

.15 

110.92 

8.57 

17.01 

.21 

699.07 

262.17 

264.06 

9 .69 

2.76 

2.58 

71.00 

3980.00 

• 14 

105.65 

8.58 

16.37 

.21 

645.25 

262.65 

264.86 

9.56 

2.54 

2.55 

72.00 

3360. DO 

.19 

100.97 

8.59 

15.79 

.21 

641.50 

263.12 

265.61 

9 .99 

2.32 

2.53 

73.00 

3290.00 

.13 

96.38 

8.59 

15.11 

.21 

637.81 

263 .56 

266.33 

9.92 

2.12 

2.50 

79.00 

3120.00 

.12 

91 .89 

8 *60 

19.99 

.20 

634.19 

264.00 

267.03 

9.35 

1.93 

2.98 

75.00 

3000. 00 

.12 

87.48 

8.61 

13.87 

.20 

630.63 

269 .92 

267.70 

9.26 

1.75 

2.95 

76.00 

2880.00 

.11 

83.16 

8.61 

13.26 

.20 

627.14 

264.82 

268.36 

9.22 

1.58 

2.93 

77.00 

2760.00 

• 11 

78.93 

8.62 

12.65 

.20 

623.70 

265.21 

268.99 

9.16 

1 .42 

2.91 

78.00 

2690.00 

.10 

79.77 

8.63 

12 . C5 

.20 

620.33 

265.59 

269.60 

9.10 

1.27 

2.38 

79.00 

2520.00 

.09 

70.69 

8.63 

11.96 

.20 

617.01 

265.96 

270.19 

9.04 

1.13 

2.36 

80 .GO 

2900.00 

.09 

66.69 

8 .64 

10.87 

.20 

613 . 7b 

266-31 

270.77 

8.99 

1.00 

2.39 

81.00 

2280.00 

.08 

62.76 

8.69 

10.28 

.20 

610.56 

266.66 

271.32 

8.93 

• 88 

2.32 

82.00 

2160.00 

.08 

58 .91 

8 .65 

9.70 

.20 

607.41 

266.99 

271.86 

8.88 

.77 

2.30 

83.00 

2090. OCi 

.07 

55.12 

8.65 

9.12 

.19 

604.32 

267-32 

272.39 

8.83 

.67 

2.28 

89.00 

1920. DO 

.07 

51 .91 

8 .66 

8.55 

.19 

601.28 

267.63 

272.89 

8.78 

.57 

2.25 

85.00 

18CD.00 

.06 

97.76 

8.66 

7.98 

.19 

598.30 

267.93 

273.38 

8 .74 

.49 

2.23 

86.00 

1680 .00 

• 0 6 

99.17 

8.67 

7.92 

.19 

595.36 

268 .23 

273.86 

8.69 

.41 

2.21 

87.00 

1560.00 

.05 

40.65 

8.67 

6.86 

.19 

592.48 

268 ,52 

279.33 

8.65 

.34 

2.19 

88.00 

1990.00 

.05 

37.19 

8 .66 

6 .31 

.19 

589.64 

268.79 

279.78 

8 .60 

• 28 

2.18 

89.00 

1320.00 

.05 

33.79 

8.68 

5.76 

.39 

586.86 

269 .06 

275.21 

8 .56 

.22 

2.16 

90.00 

1200.00 

.04 

30.45 

8 .69 

5.21 

.19 

584.11 

269.33 

275.64 

8.52 

.18 

2.19 

91.00 

1080.00 

.09 

27.17 

8.69 

9 .67 

•19 

581-42 

269.58 

276.05 

8.98 

.14 

2.12 

92.00 

960.00 

.03 

23-99 

8.69 

9.19 

.18 

578.77 

269.83 

276.45 

8.9 5 

• 10 

2 . 1C 

93.00 

890.00 

.03 

20.76 

8.70 

3.60 

• 18 

576.16 

270 .07 

276.84 

8 .91 

.07 

2.08 

99 .DO 

720.00 

.02 

17.65 

8 . 7C 

3. 08 

• 18 

573.60 

270.31 

277.22 

8 .38 

.05 

2.06 

95.00 

600.00 

.02 

19 .58 

8.71 

2.55 

.18 

571.08 

270.54 

277.59 

8.39 

.03 

2.05 

96,00 

980.00 

.02 

11.56 

8.71 

2.03 

• 1 8 

568.60 

270.76 

277.95 

8.31 

.02 

2.03 

97,00 

360.00 

.01 

8 . b 0 

8.71 

1 .52 

.18 

566.16 

270.97 

278.30 

8.28 

.01 

2.01 

98 . Od 

290.00 

.01 

5.66 

8.72 

1 .01 

.18 

563.76 

271.19 

278.64 

8.29 

• 00 

2. CO 

99.00 

120.00 

.00 

2.62 

8.72 

• 50 

• 1 8 

561.39 

271.39 

278.97 

8.21 

• 00 

1.97 

IQO.OO 

.00 

.00 

• GO 

8.72 

.00 

. 18 

559.07 

271.59 

279.29 

8 .18 

.00 


CURRENT DENSITY AT RATED TORQUE IN ROTOR BAR r 12573 . 

IN END RING = 1313R . 
IN armature = 110A7 . 


,00 


SLIP, 

PER 

CENT 


1 

1 

1 

1 

1 

1 

1 

1 

1 

20 . 1 -- 
1 
1 
1 
1 
1 
1 
1 
1 
1 

40.1 — 
1 

1 

1 

1 

1 

1 

1 

1 

1 

60. 1 — 
1 

1 

1 

1 

1 

1 

1 

1 

1 

80.1 — 
1 
1 
1 
1 
1 
1 
1 
1 
1 

100.1 — 

.00 



00 

C71 



5.00 

-1 — 
1 
1 
1 
1 
1 
1 
1 
1 
I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

-1 — 
1 
1 
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1 
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1 
1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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APPENDIX C 


PROGRAM LISTINGS 


The complete FORTRAN listings of the main program and the five subroutines, 
which together constitute the induction motor computer program, are shown in this ap- 
pendix. The main program is INDMTR and the five subroutines are, in the order given, 
CIRCT, MAGNET, SLOTS, WDGFCT, and CMBNTN. 


c 

c 


c 

c 

c 


c 


c 


c 

c 

c 

c 

c 


COMMON /CIR/ R0,R1,R2,X0,X1,X2,PW1,NSYNCH,V 1,S,I1,RPH,PF,T,HP,EFF, A 
1PIH,H1,W2,H0,FW,IMAG, V2,P0UT,PHASE A 

COMMON /MAG/ BST, BS Y , BRT , BRY, AFS T, AFSI, ATHr , ATRY, AS YD KE, A STOTH, A R Y A 
10KE, ARTOTH,LSYOKE,LRYOKE,DSS,DSa,FTOTAL,FPOLB,KSAT, AI,ATAG, ATTOr A 

COMMON /INITL/ ASTRND, AWG , CSR AP 0, DERI , DER2, D1R, D IS, D2R , D2S , D 3R , D3S A 
1,D4R,D4S,D5S,JBAR,I.B,LS,PFLUID,RTWDTH,SB,SCAREA,SFR,SPS,SKEH,SS ARE A 
2A,STNDTH,TRATED,VSCFLD,WSR1,WSR2,WSR3,iSR4, HS R5, WSS , WSS 1, NSS2, WSS3 A 
3,WSS4 A 

A 

EQUIVALENCE (RESET1 ( 1) , ASTRND) , (RESET2 (1 ) , RO) A 

A 

REAL II, IBAR, JBAR, JRTNG,LARM,LT, LTS, LTR ,LSY0KE, LR YOKE , LTOT AL, IM AG, A 
1IMAG2,NSYNCH,KPS,KDS,N,LB,NB,KRING,L,LS,KS,LREF,MATDEK,NAME A 

A 

INTEGER SSTYPE, RETYPE, ARG,RWMAT, S«MAT A 

A 

LOGICAL LAST A 

A 

DIMENSION SLIP(70), TORQUE(70), PP(61), XL3ND(3), AI(50), SMAT(13) A 
1, RHAT{13), WAREA(40), RSTVTY (5) , TMPCF(5), TITLE{13), DHSTY (5) , C A 
2LOSS(5,10), RESET1(38), RESEr2{7), NAME(13) A 


A 

NAMELIST /RATING/ NS YNCH, F,XO , X 1 , X2, RO, R1, R2, FH 1 , V 1 , Tfi ATED/STAT3R/ A 
1D,L,LTS,SPS,D0S/SSL0TS/D6S,WSS6,QS,DSS, WSS, SSTYPE ,D2S , D1S , D4S, D3S , A 
2WSS1 ,WSS2,PHIS,D5S,STWDTfl,SCAREA,CSRATO/STRWDG/CSS,PC,8,SirHAT,SPIT A 
3CH,LS, ASTRND,TSH,S,ANG,STRNDS/ROTOR/SKEW,LTR,SFR,DIR/RSLOTS/SB,D6R A 
4,WSR6,RSTYPE,D4R,D3R,WSR1 ,WSR,D2R, D1 B, HSR2, PHIR, DSR , RTHDTH/RTR WDG/ A 
5LB, NB, DER1,DER2,TER,RWMAT,TR¥,BR/AIRGAP/3,TFL0rD,VSCFLD,PFLUID, FLD A 


6NME A 

NAMELIST /FELOSS/ WCOR E, FCORE, SLOPE , BK , LT, LA ST/WNDAGE/NL, DIAREF, LR A 
1EF,RPMREF,VSCREF,C0,C1,C2,C3,C4, GAPREF, TREF, pref a 

A 

DATA RSTVTY/1,08,2.95,0.678,0.,0./,TMPCF/.00 415,.00 2, .00393,0.,0./ A 
1,DNSTY/0.0975, 0.308, 0. 32 1,0.,0./ A 

DATA (WAREA (I) ,1= 1 , 40 ) /O • 06573 , . 052 1 3 , . 0 4 134, . 0 32 78, . 02600, . 02062, A 
1.01635, .01297, .01028, 0. 008155,0.006467,0.005129,0.004 067,0.003225, A 

20.002558.0. 002028.0.001609.0.001 276.0.001012.0.0008023.0.0006363.0 A 

3.0005046.0. 0004002.0.0003173.0.0002517.0.0001996.0.0001583.0.03012 A 

455,9. 953E-5,7.894E-5,6. 260E-5,4. 964 E-5, 3. 937E-5, 3. 1 22E- 5, 2 . 476E- 5, A 
51. 964E-5, 1. 557E-5, 1.2 35E-5,9.793E-6,7.766E-6/ A 

DATA XLGND (1)/18HSLIP, PER CENT /, (PP ( I) , 1=5, 10) , KDDE/2 . , 5. , 0. , A 

120.. 4. .0.. 56/ A 

DATA BLANK/6H /, MATDEK/1 HM/ A 

A 

THE FOLLOWING ARITHMETIC STATEMENT FUNCTION GIVES THE VISCOSITY A 

OF THE FLUID IN THE MOTOR CAVITY AS A FUNCTION OF TEMPERATURE, A 
VSCSTY IS IN LBM/FT-SEC AND T IN DEG C A 

A 

VSCSTY (T)=C0 + T* (Cl + T* (C2^T* (C3 + C4*T) )) A 

A 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


10 HEAD (5, 20) DEKTYP,NAME 

20 FORMAT (A1 ,1X,13A6) 

IF (DEKTYP.EQ.MATDEK) GO TO 40 
DO 30 1=1, 13 


A 51 
A 52 
A 53 
A 54 


30 

TITLE (T)=NAME (I) 

A 

55 


GO TO 140 

A 

56 

40 

DO 50 1^1 ,13 

A 

57 

50 

SMAT(I)=HAHE(I) 

A 

58 

C 


A 

59 

c 


A 

60 

c 

IHITIALIZATION AT START OF A NEST DATA SET 

A 

61 

c 


A 

62 


LA ST=. FALSE. 

A 

63 


TBEF=25. 

A 

64 


VSCREF^O. 

A 

65 


pref=o. 

A 

66 


RPMREF=0. 

A 

67 


LREF=0- 

A 

60 


DIAREF=0. 

A 

69 


GAPREF=0. 

A 

70 


C0=0. 

A 

71 


C1=0. 

A 

72 


C2=0. 

A 

73 


C3=0. 

A 

74 


C4=0. 

A 

75 


WL=0. 

A 

76 

c 


A 

77 

c 

HEAD AND RRITE STATOR MATERIAL MAGNETIZATION DATA 

A 

78 

c 


A 

79 


READ (5,60) (AI(I) ,1=1,29) 

A 

80 

60 

FORMAT (8F10.1) 

A 

81 


WRITE (6,70) SMAT 

A 

82 

70 

FORMAT (1H1 ,4X,15HSTATOR MArERIAL/7X, 13A6) 

A 

83 


WRITE (6,80) (AI(I) ,1 = 2,29) 

A 

84 

80 

FORMAT (1HK,12X,1HB,20X,1HH//5X, 17H (KILOLINES /SQ-IN ) , 7X, 1 1 H (A-TDR N 

A 

85 


1/IN)//(F16.2,F21.2) ) 

A 

86 

c 


A 

87 


WRITE (6,90) 

A 

88 

90 

FORMAT ( 1HK,6X, 14HC0RE-L0SS D A T A /I OX , 9H CORE- LOSS , 5 X ,7 H LA M THK,5X,4 

A 

89 


1HFREQ,5X,10HFLOX DNST Y, 5X , 5H5 LOPE) 

A 

90 

C 


A 

91 

c 

READ CORE LOSS DATA FOR STATOR MATERIAL 

A 

92 

c 


A 

93 


DO 100 1=1 ,11 

A 

94 


SLOPE=0. 

A 

95 


READ (5,FELOSS) 

A 

96 


IF (LAST) GO TO 120 

A 

97 


CLOSS (1 ,1) =WCORE 

A 

98 


CL OSS (2, I) =FCORE 

A 

99 


IF (SLOPE. LT. 1 .OB-15) S LOPE= ( 1- 3 + 1 64. *LT) / ( 1 - 0 + 82. *LT ) 

A 

100 


CLOSS (3,I)=SLOPE 

A 

101 


CLOSS (4,I)=BK 

A 

102 


CLOSS(5,I)=LT 

A 

103 

100 

WRITE (6,110) WCORE,LT,FCORE, BK, SLOPE 

A 

104 

110 

FORMAT (1H ,F15. 1,F14.3,F11.1,F1 1.1,F12.1) 

A 

105 

120 

NCARDS=I-1 

A 

106 

c 


A 

107 

c 

READ ROTOR MATERIAL MAGNETIZATION DATA 

A 

108 

c 


A 

109 


READ (5,20) DEKTYP,HMAT 

A 

110 


READ (5,60) (AI( I) ,1 = 31,59) 

A 

111 


WRITE (6,130) RMAT 

A 

112 

130 

FORMAT ( 1HL,4X, 14HR0T0R M ATEBIAL/7X , 1 3 A 6) 

A 

113 


WRITE (6,80) (AI (I) ,1 = 32,59) 

A 

114 

C 


A 

115 

c 

READ WINDAGE DATA 

A 

116 

c 


A 

117 


READ (5,WNDAGE) 

A 

118 


IF (7SCREF.lt. 1.0E-15) VSCREF=VSCSTY (TREF) 

A 

1 19 

c 


A 

120 


READ (5, 20) DEKTYP, TITLE 

A 

121 

c 


A 

122 

140 

WRITE (6, 150) TITLE 

A 

123 

150 

FORMAT (1H1/1R ,2X,13A6) 

A 

124 

C 


A 

125 

c 

INITIALIZATION AT BEGINNING OF A NEW MOTOR DESIGN DECK 

A 

126 

c 


A 

127 
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DO 160 1=1 ,38 


A 128 

160 

RESET1 (I) =0. 


A 129 


DO 170 1=1 ,7 


A 130 

170 

RESET2(I)=0, 


A 131 


DSR=0. 


A 132 


DSS=0. 


A 133 


FLDNI1E = BLANK 


A 134 


PHIR=0. 


A 135 


PHIS=0. 


A 136 


RWMAT=3 


A 137 


STRNDS=1 .0 


A 138 


SW«AT=3 


A 139 


TFLOID=TREF 


A 140 


TRW=25, 


A 141 


TSW=25. 


A 142 


DO 180 1=12,61 


A 143 

180 

PP (I) =BLANK 


A 144 

C 



A 145 

C 

READ *WOTOR DESIGN* DECK 


A 146 

c 



A 147 


READ (5, RATING) 


A 148 


IP (X0+X1*X2*R0*R1*R2.GT. 1,0E'15) GO TO 720 


A 149 


READ (5, STATOR) 


A 150 


READ (5, SSLOTS) 


A 151 


READ (5,STRWDG) 


A 152 


READ (5, ROTOR) 


A 153 


READ (S,RSLOTS) 


A 154 


READ (5,PTRWDG) 


A 155 


READ (5, AIRGAP) 


A 156 

c 



A 157 

z 

RETRIEVE CORE LOSS DATA FROM ARRAY GLOSS FOR 

DESIGN LAMINATION 

A 158 

c 

THICKNESS 


A 159 

c 



A 160 


DIFF=10. 


A 161 


DO 190 1=1 ,NCARDS 


A 162 


DIFF1 = ABS(LTS-CL0SS (5,1) ) 


A 163 


IF (DIFFI.GT.DIFF) GO TO 190 


A 164 


IA=I 


A 165 


DIFF=DIFF1 


A 166 

190 

CONTINHE 


A 167 


IF (DIFF.GT,0.00 05) WRITE (6,200 ) GLOSS (5, 1 A) 


A 168 

200 

FORMAT (1HK,68HCORE-LOSS DATA IS NOT GIVEN AT 

SPECIFIED STATOR LAM 

A 169 


1INATION THICKNESS/1H ,3X,12HUSE DATA FOR,F6.3 

,12H LAMINATIONS) 

A 170 

C 



A 171 

C 

CALCULATE CORE LOSS AT DESIGN FREQUENCY 


A 172 

c 



A 173 


WFE=CLOSS ( 1, lA) * ( (F/CLOSS (2,IA) ) ♦♦CLOSS (3 ,1 A) ) 

A 174 


BK=CLOSS(4,IA) 


A 175 

c 



A 176 

c 

CALCULATE VARIOUS DIMENSIONS PROM INPUT DATA 


A 177 

c 



A 178 


DR=D-2. *G 


A 179 


T1R= (3. 1416* (DR) )/NB 


A 180 


T1S= (3. 1416*D) /QS 


A 181 


IF (SKEW,LT.1,0E-15) SKEW=AMAX1 (T1R,T1S) 


A 182 


IF (LB.LT. 1.0E-15) LB=SQRT (L*L+ SKE W* SKE W) +2.* 

(BR+TER) 

A 183 


IF (ASTRND.LT. 1 .OE-15) AST RND=W AREA (AWG) 


A 184 


SS=ASTRND*STRNDS 


A 185 


IF (SFR.LT. 1,0E-15) SFR=LTR/(LTR + 0 . 0005) 


A 186 


IF (SFS.LT.1. OE-15) SFS=LTS/ (LIS +0. 0005) 


A 187 

c 



A 188 


IF (DSS.GT.1.0E-15,OR,D1S.GT, 1. OE-15) GO TO 230 

A 189 


IF (SCAREA.GT, 1.0E-15) GO TO 240 


A 190 


IF (CSRATO.GT. 1.0E-15) GO TO 220 


A 191 


WRITE (6, 210) 


A 192 

21 0 

FORMAT ( 1HK,59HINSUFFICIENT STATOR SLOT DATA, 

SPACE FACTOR OF 0.70 

A 193 


1 ASSUMED) 


A 194 


CSRATO=0.70 


A 195 

220 

scarea=css»ss/csrato 


A 196 


GO TO 250 


A 197 

230 

SCAREA=0. 


A 198 

240 

csrato=o. 


A 199 

C 



A 200 
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250 IP (SPITCH.LE. 0.3333) KS=0. 75*S PITCH A 201 

IP (SPITCH.GT. 0.3333. AMD. SPITCH. IT. 0.6667) KS=1 .5*SPirCH-0. 25 A 202 

IF (SPITCH.GE. 0.6667) KS= 0. 75*S PITCH+0. 25 A 203 

C A 204 

CALL SLOTS ( 1 . 0, SSTTPE, WSS, HSSI , HSS2 , WSS3 ,BSS4 , HSS5 ,DSS , D1 S , D2S , D3 A 205 

1S,D4S,D5S,STHDTH,SCAREA,SSAREA,QS,D6SrHSS6, D, ICS , AXS, STWHAG, PHIS) A 206 

C A 207 

CALL SLOTS (-1 . 0 , RSTYPE, WS E, HS R 1 , WS R2, WSH 3, WS S4, 5TSR 5, DSR, D 1R ,D2H , D A 208 

13H,D4R,0.,RTHDTH,SB,RSAREA,NB,D6R,HSR6rDRr1 . 0 , AXR , RTHHAG, PHIR) A 209 

C A 210 

C STATOR AND ROTOR IRON WEIGHTS A 211 

C A 212 

HSTOTH= (3.1416* (D+DSS)*DSS- (SSAREA+QS) ) *L*SPS*. 283 A 213 

HSYOKE= (0.7854*(DOS*DOS-(D*2.*DSS) **2) ) *L*SFS*0.283 A 214 

HSTAT=RSTOTH+HSYOKE A 215 

HROT= (0.7854* (DR *DR- DIR *DIR) -NB* (RSAHBA) ) *L*SPR*0.2B3 A 216 

C A 217 

C END RING DIMENSIONS A 218 

C A 219 

IF (DER1.LT. 1.0E-15) DER 1 =DR- 2 . * (D4 R + D3R) A 220 

IF (DER2.LT. 1.0E-15) DER 2 = AM A XI ( DR-3 , *DSR ,1 . 1 *DI R) A 221 

SER=0.5* (DER1-DER2)*TER A 222 

C A 223 

IF (CSRATO.lt. 1.0E-15) CS RATO=CS5 *SS/SC AR EA A 224 

P=FLOAT(IFIX(((120.*F)/NSYNCH)*0.1)) A 225 

N=(QS*CSS)/(PC*3.) a 226 

DBS=(DOS-D)*0.5-DSS A 227 

DBRS= (DR-DIR) *0. 5-DSR A 223 

C A 229 

C CHECK IF STATOR-ROTOR SLOT COMBINATION IS ACCEPTABLE A 230 

C A 231 

CALL CMBNTN (QS,NB,P) A 232 

C A 233 

C CALCULATE DISTRIBUTION AND PITCH FACTORS A 234 

C A 235 

CALL WDGFCT (60 . , P , Q S, KDS , PC , KP S , SPI TCH ) A 236 

C A 237 

C CARTER COEFFICIENTS AND EFFECTIVE AIRGAP A 238 

C A 239 

IF (RSTYPE.GT. 2) GO TO 260 A 240 

CCR=(T1R* (5.*G+WSR))/(T1R*(5.*G+WSR) -WSR*HSR) A 241 

GO TO 270 A 242 

26 0 CCB=(T1R*(4.4*G^'0.75*HSR1) )/(TlR*(4. 4*G+0 . 75*WSR 1 ) - WSR 1 *W5R 1 ) A 24 3 

270 IF (SSTYPE.GT. 2) GO TO 280 A 244 

CCS= (T1S*(5. *G4-WSS) ) /(T1S*(5. *1<^ WSS) -HSS*WSS) A 245 

GO TO 290 A 246 

280 CCS=(T1S*(4.4*G + 0.7S*WSS1) )/(TlS* (4.4*G + 0.75*WSS1) -WSS1*HSS1) A 247 

290 GE=G*CCR*CCS A 248 

C A 249 

C STATOR RESISTANCE CALCULATION (R1) A 250 

C A 251 

IF (SSTYPE.BQ.2.0R.SSTYPE.EQ.4. OR.SSTYPE. EQ- 6) S ALP HA = { 0 . 5* ( WSS 4^^ W A 2 52 

1SS5) +S-2.*HSS6)/(3. 1416*(D+DSS) /QS) A 253 

IF (SSTYPE.EQ.1.0R.SSTYPE.EQ.3) S ALPH A= (»SS +S-2. *WSS6 ) /T IS A 254 

IF (SSTYPE.EQ.5) SALPHA = ( WSS3 + S- 2 . ♦ WSS6 ) / (3 . 14 16 * ( 2. *D4S ♦ HSS 3) /Q A 255 

IS) A 256 

CALPHA=SQRT (1.-SALPHA**2) A 257 

AY=(3. 1416*(D+DSS) *SPITCH)/(P*CALPHA) A 258 

IF (LS. LT. 1 .OE-15) LS= A Y+ 2 . *B+D5 S +L A 259 

IF (R1. GT. 1. OE- 15) GO TO 300 A 260 

R1=(LS*N*RSTVTY(SWMAT) *1. OE-6) /(PC*SS) A 261 

R1=R1* (1 .+THPCF (SWMAT)* (TSH-20. ) ) A 262 

C A 263 

C AXIAL EXTENSION OF END TURN AND OVERALL ARMATURE LENGTH A 264 

C A 265 

300 ENDTRN=^A Y*0.5*S ALPHA+B+DSS A 266 

LT0TAL=L+2.*ENDTRN A 267 

C A 268 

C ARMATURE WEIGHT AND TOTAL WIRE LENGTH A 269 

C A 270 

LARM= (LS*CSS*QS*STRNDS) /12. A 271 

WARM=DNSTY(SWMAT) *LARM*ASTRND*12 . A 272 

C A 273 
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C ROTOR RESISTANCE CALCDLATION (R2) A 274 

C A 275 

RRSTVY=1.0E-6*RSTVTY (RWMAT) ♦(1. 0+TMPCF (RWMAT) ♦ (TRH-20. ) ) A 276 

IF (R2.GT, 1 .OE-15) GO TO 310 A 277 

RATIO=DER2/DER 1 A 278 

KRING = 0,50*P* (1. -RATIO) * ( 1 • +R AT IO**P) / ( 1 . -R ATIO**P) A 279 

AY=( (N*KPS*KDS)**2)»3.*PRSTVY A 280 

R2BAR=AY*((LB-TER)/(SB*NB) ) A 281 

R2RING=AY* ( (0. 64*DBR 1 * KRI NG) / {?♦ P*S ER) ) A 282 

R2=R2BAR+R2RING A 283 

C A 284 

C ROTOR WINDING WEIGHT AND TOTAL ELECTROMAGNETIC MOTOR WEIGHT A 285 

C A 2 86 

31 0 HRWNDG=DNSTY (RWMAT) ♦ (NB^SB* (LB- 2 .♦TER) + SE P* 3 . 1 4 1 6^ { DER 1 ♦DER2) ) A 287 

WEIGHT=WARM+ WRWNDG+WROT+WSTAT A 283 

C A 289 

C MAGNETIZING REACTANCE (AIR GAP ONLY) A 290 

C A 291 

X0AG = 7. 66E-7*F*(N/2. *KPS*KDS) *D ♦!/ (P *P*GE) A 292 

C A 293 

C PRIMARY SLOT LEAKAGE REACTANCE (XSS) A 294 

C A 295 

AY=6. E-7^N*N*F*L A 296 

XSS=AY*AXS/QS A 297 

C A 298 

C SECONDARY SLOT LEAKAGE REACTANCE (XRS) A 299 

C A 300 

XRS=AY^ (KPS^KDS) ♦*2*AXR/NB A 301 

C A 302 

C ROTOR AND STATOR END-CONNECTION LEAKAGE REACTANCE A 303 

C A 304 

DC= (D + 2.^ (D4S+D3S+D2S) f DIS-n. 50 ♦ (DEKUDER2) ) *0. SO A 305 

AY:= ( (N^N*F^ (KPS*KDS) ♦♦2) /P) ♦2.4E-7 A 306 

F 1 = 1 . 570 8* (D+DSS) ♦$ PITCH ♦SQRT (1 . 0-C ALPH A*CA LPHA ) / (P*C ALPHA) A 307 

XSE = AY^ (B+0.5* (F1+DSS/2. ) ) A 308 

XPE= (0. 72 5 + AY/P) * ( 2 . *P ♦BR + ( 3 . 1 4 1 6*D*DC ) / ( 1 . 7^TE R+ 0 . 6 * (DER1-DER2) +1 A 309 
1 .4*DC) ) A 310 

C A 311 

C SKEW REACTANCE (XSK) A 312 

C A 313 

XSK = 0. 5* (XOAG/12.) ♦ (P*SKEW/D) * + 2 A 314 

C A 315 

C STATOR AND ROTOR ZIGZAG LEAKAGE REACTANCE A 316 

C A 317 

XZ=0. 8333*X0AG+ (KS/ (KP S*KDS) **2) A 318 

XPZ=XZ*((6./CCR-1.)/(S.^(NB/P)**2)) A 319 

XSZ=XZ^ ( (6./CCS-1. )/ (S. ♦ (QS/P)**2) ) A 320 

- A 321 

C PERIPHERAL AIR-GAP LEAKAGE REACTANCE (XP) A 322 

C A 323 

XP = 0.52S+X0AG* (P*G/D) * + 2 A 324 

C A 325 

C TOTAL ARMATURE AND ROTOR LEAKAGE REACTANCES (XI AND X2) A 326 

C A 327 

IF (ABS (X 1) .LT. 1. OE-15) X 1 = XSS+ X SE + XSK+ XSZ^ XP A 328 

IF (X2. LT. 1 .OE-15) X2=XRS +XRE+X5 K+XKZ A 329 

C A 330 

C WRITE OUTPUT A 331 

C A 332 

WRITE (6,320) NS Y NCH , F , P , V 1 A 333 

320 FORMAT ( 1 HL , 5X , 6HR ATING/1 OX, 1 7H S Y NCH RONOUS SPEED ,F2 6 . 0 , 4 H RPM/10X, A 334 

19HFRBQUENCY,P34.0,3H HZ/ 1 OX , 5 HPO L ES , F38 . 0/ 10X, 1 1 HL- N VOLTA GE, F33 . 1 A 335 

2,5H VOLT) A 336 

IF (TRATED.GT. 1 .OE-15) WRITE (6,330) TRATED A 337 

330 FORMAT (1H ,9X,12HRATED TORQUE, F 32 . 1 , 7H IN-LBS) A 338 

C A 339 

WRITE (6,340) D, DOS, DB S, L ,LTS , SF S, WSTA T A 340 

340 FORMAT ( IHL, 5X, 6HST ATOR/1 OX, 1 3HBORS DIA METE R , F33 . 3/1 0 X , 1 6 HOUTSIDE A 341 

1DIAMETER,F30.3/10X,16HDEPTH BELOW S LOT, F30. 3/ 10 X, 6HLEN GTH, P 40. 3/1 0 A 342 
2X, 20HLAMTNATION THICKNESS , F26 . 3/1 OX , 1 5H ST ACKI NG F ACTOR , F3 1 . 3/ ID X , 1 A 343 

38HSTATOR IRON WEIGHT , F28 . 3) A 344 

C A 345 
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WRITE (6,350) SSTYPE,QS 

350 FORMAT ( 1HL, 5X, 1 2HSTATOR SLOTS/1 OX, 9BSLOT TYPE,I8 ,9X, 1 2HNO. OP SLO 
1TS,F5,0) 

IF (SSTYPE. EQ. 1.0R. SSTYPE.EQ.3) WRITE (6,360) HSS,DSS 
IF ( (SSTYPE/2) ♦2,EQ, SSTYPE) WRITE (6,370) STWDTH, DSS 
IF (SSTYPE. EQ. 5) WRITE (6,380) WSS3, DSS 
360 FORMAT (10X,10HSLOT WIDTB,F10. 3, 6X, 10HSLOT DEPTH,F10.3) 

370 FORMAT (10X,11HT00TH WIDTB,F9 . 3, 6X, 1 OHSLOT DEPTH,F10.3) 

380 FORMAT (10X,13HSLOT DI A METER, F7. 3, 6X , 10 HSLOT DEPTH,F10.3) 

WRITE (6,390) WSS 1, D 1 S, WSS2, D2S, HSS3 , D3 S, WSS4 ,D4 S , WSS 5, D5S , HSS6, D6 
1 S,SCAREA,SSAREA,CSRATO 

39 0 FORMAT (1 OX ,4 HWSS 1 , FI 6 . 3 ,6X, 3HD 1 S , F 17. 3/ 10X, 4HWSS 2, F 1 6 . 3, 6X, 3HD 2S, 
1F17. 3/1 OX, 4HWSS3, F16. 3,6X,3HD3S,F17.3/10X,4HHSS4,F16. 3 ,6X,3HD4S,Fl 

2 7.3/10X,4HWSS5,Fl6. 3 ,6X,3HD55,F1 7. 3/1 OX, 4 HWSS 6, FI 6. 3, 6X, 3RD6S, F 1 7. 
33/10X, 11HUSABLE AR EA , F9 . 3 , 6X, 1 OH TOTA L A RE A, F 1 0 . 3 /I 0 X, 1 2HS P ACE FACT 
40R,F8.3) 

C 

WRITE (6,400) SWMAT,CSS,PC,SPITCH, B,SS, ASTRNO, LS, S, TSW, ENDTRN, LTOT 
1AL,KPS,KDS, WARM, L ARM, STRNDS 

400 FORMAT (1HL,5X, 14HST ATOR WINDING/ 1 0X, 8 HMATER lAL, 135/1 OX, 19HC0ND0CT 
10RS PER SLOT,F24.0/10X, 17HPARALLEL Cl RS UI TS ,F26 . 0/1 OX , 5 HPITCH, F4 1 . 
23/10X,27HAXTAL EXTENSION BEYOND CORE, Fl 9. 3/ 10X, 2 3HCON DOCTOR CROSS- 
33 ECT ION, E27. 3/10X,20HSTRAND CROS S- SBCTI ON ,E30 . 3 /I 0 X , 1 6 HCONDOCTOR L 
4ENGTH,F30.3/10X,21HCLRNCE BTWN E ND-TURNS , F25. 3/ 10X, 15HTEMPER ATORE 
5 (C) , F28. 0/10X, 21HAXIAL END-TURN LENG TH , P2 5. 3/1 0 X , 23 R07ER ALL APMATU 
6RE LENGTH, F23.3/10X,12HPITCH F ACTOR, F34 . 3/ 1 OX, 19 HDISTR IBOTION FACT 
70R, F27, 3/10X, 15HARM ATORE WEIGHT, F3 1 . 3/1 OX ,2 6 RTOTAL ARMATURE WIRE L 
8EN3TH,F20.3 ,5H FEET/1 0 X , 1 7 HST R A N DS/CON DUCTO R, F26 . 0) 

IF (AWG. GT.O) WRITE (6,410) A WG 
410 FORMAT ( 1H ,9 X, 1 1HSTR A ND SIZE, 132) 

C 

WRITE (6, 420) DR, DIR , LTR , SFR , SK E W , DB RS , WROT 
420 FORMAT (1 H 1 , 5 X ,5R ROTOR/ 10 X , 1 4 HBOT OR DI A H ET ES, F3 2. 3/ 10 X , 1 5H IN SIDE D 
1 lAMETER, F31. 3/10X, 20HI.A M INA TI ON TH TCKNB SS , P26 . 3 /I 0 X , 1 5 HST ACKTNG FA 
2CT0R,F31 ,3/10X,9HSLOT S KEW , F3 7. 3 / 1 0X , 1 6 HD EPT R BELOW S LOT, F 30 . 3/1 OX 
3, 17HR0T0R IRON H EIGHT , F 2 9 . 3) 

C 

WRITE (6,430) RSTYPE,NB 

430 format (1HL,5X,11HR0T0R S LOTS/1 0 X , 9 HSLOT TY P E, I 8, 9X , 1 2HNO , OF SLOT 
IS, F5.0) 

IF (RSTYPE.EQ. 1- OR. RSTYPE. EQ. 3) WRITE (6,440) WSR,DSR 
IF ( (RSTYPE/2) *2. EQ. RSTYPE) WRITE (6,450) RTHDTH,DSR 
IF (RSTYPE.EQ. 5) WRITE (6,460) MSR3,DSR 
440 FORMAT (10X,10HSLOT W TDT H , F 1 0 , 3 , 6X , 1 OHS LOT DEPTH,F10.3) 

450 FORMAT ( 10X, 1 1 HTOOTH WIDTH, F9 . 3 , 6 X, 1 Ofl SLOT DEPTH , FI 0 , 3 ) 

460 FORMAT (10X,13HSLOT DT A METER, F7. 3, 6X , 1 0 HSLOT DEPTH,F10.3) 

WRITE (6, 470) W SR 1 , D 1 R , WSR 2, D 2R , W SR 3 , D3 R , WS R4 , d4 R , WSR 5, D6P, WSR6, SB 
1 ,RSAREA 

47 0 FORMAT (10X,4RWSR 1 , F 1 6 . 3 , 6X , 3 HD 1 R, F 1 7. 3/ 1 0X , 4HWSR2, FI 6. 3, 6X, 3RD2R, 

1F17, 3/10X, 4HWSR 3, F 1 6 . 3 , 6 X , 3HD 3R , F 1 7 . 3 /1 0 X ,4H WSR4 , F1 6 . 3 ,6X,3HD4R,F1 
27.3/10X,4H WSR5,F16. 3,6X,3HD6R,F1 7. 3/ 1 0X , 4 HW S R6, F 1 6. 3/ 1 0X, 1 IHUSABLE 

3 AREA, F9. 3, 6X, 10HTOTAL AREA,F10.3) 

C 

WRITE (6, 480) R WM AT, LB , SB , DER 1 , D ER2 , TER , B R , TRW , WR WNDG , R2 B A R, R2 RING 

48 0 FORMAT ( 1 HL ,5 X , 1 3HROTOR W I NDI NG/ 1 OX , 8 HM AT ERIAL, 135/1 OX, 10HBAR LENG 

1TH,F36. 3/10X,17HBAR CROSS-SECTION, F29. 3/ 10X, 20HEND-RING OUTSIDE DI 
2A,F26.3/10X, 19HEND-RING INSIDE D I A , F27 . 3/10 X, I 8HE ND-RI NG THICKNESS 
3,F28.3/10X,24HSTACK-TO-END-RING CLRNCE, F22. 3/10X, 23HWINDING TEMPER 
4ATURE (C) ,F20.0/10X,6HWEIGHT,F40.3/10X,27HCOMPONENT OF R2 DUE TO B 
5ARS,F19.3/10X,32HCOMPONENT OF R2 DOE TO END RINGS,F14.3) 

C 

WRITE (6,490) G,GE,X0 AG 

49 0 FORMAT (1 HL , 5X, 6 HAIR GAP/1 0 X, 1 3H ACT UAL AIRGAP, F34. 4/10X, 16REFFECTIV 

IE AIRGAP, F31. 4/10X, 36HMAGNETIZING REACTANCE (AIR GAP ONLY),F9.2) 

C 

WRITE (6, 500) XSS,XRS, XSE , XRE , XS K , XSK , XSZ , XR Z, XP 
500 FORMAT (1HL,5 X, 24 H LEAKAGE REACTANCES (OHM) / 33X, 6HSTAT0R , 1 IX, 5RROTO 
1B/10X, 4HSLOT,F26. 3, FI 6. 3/1 OX, 1 4H END-CON NECTI ON , F1 6. 3,F16.3/10X,4HS 
2KEW,F26,3,F16.3/10X,7HZIG-ZAG,F23. 3, FI 6. 3/1 OX, 10HPER IPHEB A L, F20. 3) 
C 

WRITE (6,510) WEIGHT 

510 FORMAT (lHL,5X,6HWETGHT/10X,23HrOTAL (ELECTROMAGNETIC) ,F23.3/1H1) 

C 
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C CPOSS-SECTIONAL AREAS AND LENGTHS OF FLUX PATHS NEEDED FOR A 419 

C HAGNETIC CALCDLATIONS A 420 

C A 421 

ASYOKE=DBS*L*SFS A 422 

LSYOKE=3.1416* (DOS+D+2.*DSS)/(4.0»P) A 423 

ARTOKE=DBRS*L*SFR A 424 

LR yOKE=3 . 1 4 1 6 * ( DR-2 . ^DS R-^ DIR) / ( 4 . 0 *P) A 4 25 

ARTOTH=RTHMAG*L»SFR*NB A 426 

ASTOTH=STHnAG*L*SFS*QS A 427 

C A 428 

C NO-LOAD MAGNETIC CALCDLATIONS A 429 

C A 430 

XX=1,0 A 431 

XY=1.0 A 432 

IF (XO.GT. 1.0E-15) XX=0.0 A 433 

IF (R0,GT.1.0E-15) XY=0.0 A 434 

X0=X0+(0.5*X0AG) *XX A 435 

W0={WSYOKE+HSTOTH) *HFE*3.0 A 436 

R0=(5.*V1*V1/W0)*XY+R0 A 437 

S=0. A 438 

ICNT2=0 A 439 

520 ICNT2=ICNT2+1 A 440 

IF (ICNT2.GE. 16) GO TO 550 A 441 

ICNT1=0 A 442 

530 ICNT1=ICNT1+1 A 443 

IF (ICNTl.GE.il) GO TO 540 A 444 

call CIRCT a 445 

R0OLD=R0 A 446 

FTOTAL=V2*P*1.0E4'05/(1.414+N*F*KPS*KDS) A 447 

FPOLE=FTOTAL*0. 637/P A 448 

BG~FT0TAL/(3. 1416*D*L) A 449 

ATAG=BG*GE*313. a 450 

CALL MAGNET A 451 

H0= (WSTOKE* (BSY/BK)**2 + HSTOTH* (BST/BK) **2) *WFE*3.0 A 452 

RO=((3.*V2*V2/WO)-RO) *XY+R0 A 453 

IF (ABS (R0-R0OLD)/R0. GE. 0-001) GO TO 530 A 454 

540 IMAG2=2.22*P»ATT0T/(3,*N*KPS*KDS) A 455 

X0^X0+ ( (V2/ (0.5* (IMAG + IMAG2) ) ) -X0> *XX A 456 

IF (ABS ( (IMAG-IMAG2) /IMAG) -GT.. 005) GO TO 520 A 457 

550 CURDEN=(SQRT(IHAG**2+(V2/R0) **2) )/(PC*SS) A 458 

IF (ICNT1.GE.11) WRITE (6,560) A 459 

IP (ICNT2.GE. 16) WRITE (6,570) A 460 

IF (KSAT.EQ.O) WRITE (6,580) A 461 

560 FORMAT (1H ,38HSHUNT RESISTANCE RO FAILED TO CONVERGE//) A 462 

570 FORMAT ( 1H , 38HM AGNETIZI NG CURRENT FAILED TO CONVERGE//) A 463 

580 FORMAT (1H ,17HWACHINE SATURATED//) A 464 

C A 465 

C WRITE RESULTS OF NO-LOAD MAGNETIC CALCULATIONS A 466 

C A 467 

WRITE (6,590) SMAT, AI (1 ) , BK, W FE A 468 

590 FORMAT (1H , 5X , 1 7HSTA TOR MATERIAL -,1R , 1 3A6/24X,7HB MAX =,F5.0/24 A 469 

1X,12HCORE LOSS AT,P6.1,10H KL/S Q-IN=, F5. 1 , 5H W/LB) A 470 

WRITE (6,600) RHAT,AT(31) A 471 

600 FORMAT ( 1H K ,5X , 1 7H ROTOR MATERIAL — ,1H , 1 3A6/24 X , 7 HB MAX =,F5.0) A 472 

C A 473 

C WRITE NO-LOAD MAGNETIZATION CHARACTERISTICS A 474 

C A 4 75 

WRITE (6,610) FTOTAL,FPOLE,BG,BST,BSY,BfiT,BRY,ATAG,ATST,ATSY, ATRT, A 476 

1 ATRY,ATTOT,TMAG,V2,CURDEN,W0 A 477 

610 FORMAT ( IHL, 5X, 29HMAGNETIZATION CHARACTERISTICS/7 X25H (NO-LOAD, RA A 478 

1TED V0LTAGE)//9X,18H TOTAL USEFUL FLUX, F28. 2, 1 OH KILOL INES/9 XI 7 H U A 479 

2SEFUL FLUX/P0LE,F29.2//9X15H FLUX DE NSITI ES/ 1 3X7H AIR G AP, F35. 2, 9H A 480 

3KL/SQ-IN/13X 13H STATOR TOOTH, F29 . 2/1 3 XI 2H STATOR YOKE , F30 . 2/ 1 3X 1 2H A 481 

4 ROTOR TOOTH, F30.2/13X11H ROTOR YOKE, F3 1 . 2//9X 22H AMPERE-TURNS PER A 482 

5 POLE/13X7H AIR GAP , F3 5. 2, /I 3X 1 3H STATOR TOOTH ,F29 . 2/1 3 XI 2 H STATOR A 483 

6YOKE,F30.2/13X12H ROTOR T OOTH, F3 0 • 2/ 1 3X 1 1 H ROTOR YOKE , F 3 1 . 2//1 3X6H A 484 

7 TOTAL, F36-2//9X20H MAGNETIZING C URRENT,F26 .2 ,8H AMPERES/10X, 14HAI A 485 

8RGAP VOLTAGE,F31-2/10X,20HN-L. CURRENT DENS XTY, F25. 2/ 1 0X, 9HCORE LO A 486 

9SS,F34.0,5H WATT) A 487 

C A 488 

C SCALE WINDAGE LOSS FROM REFERENCE CONDITIONS TO DESIGN CONDITIONS A 489 

C A 490 
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WRITE (6, 620) 

620 FORMAT (1 H L ,5 X ,7 H WI ND AGE) 

IF (FW1.GT.1.0E-15) GO TO 700 
IF (WL.LE. 1 .OE-15) GO TO 630 

IF (DIAREF*LREF*RPMREF*GAPREF.GT. 1.0E-1 5) GO TO 650 
630 WRITE (6,640) 

640 FORMAT (1 HK , 39 HIN SO FFIC lENT DATA TO SCALE WINDAGE LOSS//) 

GO TO 700 

65 0 FW1 = WL* ( (DR/DTAREF) **3. 25)* (L/LREF) ♦ ( ( NSY NCH/RPMR EF) ♦♦2. 5) ♦ ( (GAPRE 

1F/G) ♦♦0.25) 

IF (VSCREF.LT. 1.0E-15) GO TO 670 
IF (VSCFLD.GT. 1. OE-15) GO TO 660 
VSCFLD = VSCSTY (TFLUID) 

IF (7SCFLD.LT. 1. OE-15) GO TO 670 
660 FW1=FW1^ { (VSCFLD/VSCREF) ♦♦0.50) 

670 IF (PREF.LT. 1. OE- 15) GO TO 680 

IF (PFLUID.lt. 1. OE-15) GO TO 680 
FH1=FW1+ (PFLDTD/PRFF) 

C 

C WRITE WINDAGE DATA 

C 

680 WRITE (6,690) FW1 , WL, DR, DI AREF, L, LREF, NSYNCH, RPMREF,G, GAPREF,TFLOI 
1 D,TREF, VSCFLD, VSCR EF, P FLUID, PREP , FLD NME 
690 FORMAT (1H , 3 1 X , 6HDES IG N, 9X, 9 HRE FERENCE/26X , 2 ( 6X , 9HC0N DITTON ) //I OX 
1, 15HWINDAGE LOSS, W, F13. 0, F1 5 . 0/1 OX , 8HDI A METER , F23 . 3 , F 1 5. 3/10X, 6HL 
2ENGTH,F25.3,F15. 3/10 X, 3H R PK, F25 , 0 , F1 5 . 0/1 OX , 3 HG AP , F23 . 3 , F 1 5. 3/1 OX, 
311HTEMP, DEG C,F17.0,F15.0/10X, 18HVSCSTY, LBH/FT-SEC, E 1 3. 3, E 15. 3/1 
40X, 19HPRESSURE, LBS/SQ- IN , F1 1 . 3 , F 1 5 . 3//1 0 X, 5 HF LUI D , 26 X , A6 /1 H1) 

GO TO 720 

700 WRITE (6,710) FW 1 

710 FORMAT (1H ,9 X , 3 3 H WI NO AGE LOSS AT SYNCHRONOUS S P SED, F 1 0 . 0 , 5H WATT/ 
1 1H1) 

C 

C WRITE VALUES OF EQUIVALENT CIRCUIT ELEMENTS 

C 

720 WRITE (6,730) P 1 , X 1 , R 2, X2 , R 0, XO 

730 FORMAT (1 H K ,5 X , 29 HEQUI V ALENT CIRCUIT PA R AM ET ERS/ 10X , 4 HR 1 =,F9.3,15 
1X,5HX1 = , F7. 3/10X, 4HR2 = ,F9. 3, 1 5X,5HX2 = , F7 . 3 /I 0 X ,4 H RO =,F9.3,15 
2X,4HXO =,Ffi,3) 

C 

C EQUIVALENT CIRCUIT ANALYSIS 

C 

KT=1 

IF (TRATED.LT. 1. OE- 15) KT=3 

DELTAS^I -0 

SMAX=100. 

IA=TFIX ( (SMAX/ (50. + DELTAS) ) +0.5) 

S = 0. 

1 = 0 

TOLD=0- 

WRITE (6,740) VI ,F 

740 FORMAT ( 1 HK , 5X , 2 OHMO TOR PERFORMANCE AT,F7-2,6H VOLT , , F7. 1 , 3 R HZ//6 
1 X,6HT0RQUE,4 X,4HSLTP,6X,3HRPM, 12X,5HP-0UT, 12X, 1HI,7X, 3 HEFF , 7X , 2H PF 
2, 6X, 4HP-IN, 6X, 3HPRI, 6X, 3HSEC,7X, 4HIRON, 9X ,2HFW/5X,8H ( I N-LBS) , 1X,9H 
3 (PERCENT) , 14X,4H (HP) ,4X,6H (WATT) ,5X, 5H (AMP) , 2X, 9H (PERCENT) , 1 1H (*=L 
4EADING) ,6H (WATT) , 5 X , 4 HLOSS , 5 X , 4 H LOSS , 6X , 4 HL OSS , 7 X , 6R ( W ATT) /9 3X,6R ( 
5WATT) , 3X, 6H (WATT) ,4X,6H (WATT) //) 

C 

750 S=S+DELTAS 
1=1+1 

IF (S.GT.SMAX) GO TO 870 
760 CALL CIRCT 

IF (T.GT. 1 .OE-15) GO TO 780 
WRITE (6,770) S 

770 FORMAT ( 1H ,5X,44HF+W TORQUE EXCEEDS AVAILABLE SHAFT TORQUE AT,F8. 
13,13H PERCENT SLIP) 

IF (S.GT. 15.) GO TO 870 
T=0. 

GO TO 820 

780 GO TO (790,840,000) ,KT 

790 IF (T. GE.TRATED) GO TO 830 

TOLD=T 
SOLD=S 
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80 0 WRITE (6,810) T, S, RPM, H P, FOOT , II , BFF, PF, PH AS E, P IN , W 1, W 2, WO , FW A 564 

810 FORMAT (1H , 7 F 1 0 . 2, F8. 2, A 1 , F 9. 2 , 4F 1 0 .2) A 565 

820 J=I/IA A 566 

IF (J*IA.NE,I) GO TO 750 A 567 

IF (J.GT.50) GO TO 750 A 568 

SLIP(J)=S A 569 

TORQUE (J)=T A 570 

GO TO 750 A 571 

C A 572 

C CALCULATE VALUE OF S AT TORQUE TRATED A 573 

C A 574 

830 S={ ( (TRATED-TOLD) /(T-TOLD) )* (S-SOLD) ) +SOLD A 575 

KT=2 A 576 

GO TO 760 A 577 

C A 578 

C WRITE MOTOR CHARACTERISTICS AT RATED TORQUE A 579 

C A 580 

840 IF { (ABS (T-TF ATED) ) . GT.0. 005) GO TO BJO A 581 

WRITE (6,850) A 582 

850 FORMAT (1H ) A 583 

WRITE (6,810) T, S,RPM, HP,POUT,I 1 ,£FF,PF,PH ASE,PIN, W1 ,H2 , WO , FW A 584 

WRITE (6,850) A 585 

C A 586 

C CURRENT DENSITIES AT RATED TORQU F. A 5 87 

C A 588 

IF (NB.LT.1.0E-15) GO TO 860 A 589 

WBAR= (W2/NB) * (R2BAP/R2) A 590 

IBAR = SQRT{ (WBAR + SB) /(RRSTVY + LP) ) A 591 

JBAR=IBAR/SB A ^92 

WRING= (W2-WBAR) *0.5 A 593 

JRING = SQRT (WRI NG/ (R RS T V Y* S ER* 1 . 5 70R* (DEK1 + r)KR2) ) ) A *^94 

CURDEN=I1/(PC*SS) A 595 

860 KT=3 A 596 

S = SOLD + DEL'^AS A 597 

GO TO 760 A 598 

C A 599 

870 IF (JBAR.GT.1 .OE^15) WRITE (6,880) J BAR , J R I*; G, CHR DCK A 600 

88 0 FORMAT ( 1H K, 1 0X , 46H CUP R EN T DENSITY AT RATED TORQUE IN ROTOR BA'^ A 601 

1F7.0/44X,13HIN END PING = , F7 . 0/4 4 X , 1 3 HI N ARMArURE =,F7.0) A 602 

C A 60 3 

C PLOT TORQUE-SPEED CURVE A 604 

C A 605 

PP(1)=J A 606 

IA=J/2+10 A 607 

PP (I A) =XLGND (1 ) A 608 

PP(IA+1) =XLGND(2) A 609 

PP (IA + 2) =XLGND (3) A 610 

CALL PT.OTXY ( SL TP , TO PQ UE , KODE , ? P ) A 6 11 

WRITE (6,890) A 612 

890 FORMAT ( 2HPL , 70X , 1 4HTnRQ UE , IN-LBS) A 613 

GO TO 10 A 614 

END A 615- 


C 


SUBROUTINE CIRCT 

COMMON /CIR/ R0,R1,R2,XD,X1,X2,FW1 , NSYNCH , V 1 , S , I 1 , R PM , PF,T , HP, EFF, 
1PIN,W1 ,W2,H0,FW,IMAG,V2, POUT, PHASE 
REAL NSYNCH, II, 12, IMAG 
COMPLEX D, ZO , Z1 ,Z?,E1 ,E2, IA,IB, IC 
DATA STAR,BLANK/1H+, 1H / 

C=2, 5 

PHASE=BLANK 

POUT=0. 

EFF=0- 
HP = 0. 

T=0. 

E1=CMPLX (VI, 0.) 

Z1=CMPLX (R1 ,X1 ) 

Z0 = CMPLX (RO*XO*XO/(RO*RO + XO*XO) , X 0*R 0* R 0/ (RO *R0+ XO* XO ) ) 


P 1 
B 2 
B 3 
B 4 
B 5 
B 6 
B 7 
B 8 
P 9 
B 10 
P 11 
B 12 
R 13 
R 14 
B 16 
B 16 




IF (S.LT.1.0E-10) 00 TO 10 

B 

17 


Z2=CMPLX (R2*1 OO./S, X2) 

R 

18 


D= (Z 1 + ZO) * (Z2+Z0) -ZO*ZO 

B 

19 


IA=(E1* (Z04-Z2) )/D 

B 

20 


IB= (E1*Z0) /D 

B 

21 


IC=IA-IB 

B 

22 


GO TO 20 

B 

23 

10 

IA=E1/(Z1+Z0) 

B 

24 


IB= (0. , 0- ) 

B 

25 


IC = IA 

B 

26 

20 

E2= (lA-IB) *Z0 

B 

27 


A=RRAL (E2) 

B 

28 


B = ArMAG (E2) 

B 

29 


V2=SQRT { A*A+B*3) 

B 

30 


IMAG=V2/X0 

B 

31 


WO=V2*V2/R0^3. 

B 

32 


A = RS AL (lA) 

B 

33 


B= AIM AG (lA) 

B 

34 


IF (B.GT.O.) PHASE='STAP 

B 

35 


I1 = S0RT (A*A+B*B) 

B 

36 


PF=^A/I1 

B 

37 


AUREAL (IB) 

B 

38 


B^AIMAG(IB) 

B 

39 


I2 = SQRT (A* A+R + B) 

B 

40 


H 1 = I1*n*R 1*3. 

B 

41 


W2=I2*T2*R2*3. 

B 

42 


RPM^NSYNCH* (1. -S/100.) 

B 

43 


FW=FW1* (FPM/NS YNCH) **C 

B 

44 


PIN= VI *r 1*PF*3 . 

R 

45 


IF (S.GT.O.) GO TO 30 

B 

46 


POUT=-F W 

B 

47 


GO TO 40 

B 

48 

30 

POUT= W2* ( (1 00 . - S) /S) -FW 

B 

49 

40 

EFF=100. *POOT/PrN 

B 

SO 


HP=POUT/746. 

B 

51 


IF (S.LT.99.9) GO TO 59 

B 

52 


T= { (. 8 47E4/S) * W2) /N SYNCH 

B 

53 


GO TO 60 

R 

54 

50 

T= (HP/PPM) *6. 34E4 

B 

55 

60 

RETURN 

B 

56 


END 

B 

57- 



SURROUTINE MAGNET 

C 

1 


COMMON /MAG/ BST, RS Y, BRT , 0 RY , AT S P , AT S Y , ATRI, ATRY, AS Y3 K E, A STOTH, A R Y 

c 

2 


10KE, ARTOTH^LSYOK E,LR YOKF, DSS,DSR , F TOTA L , F PO LE , KS AT , AT , ATAG, A'^’TOP 

c 

3 


DIME NSION AI (60) 

c 

4 


REAL LSYOKE,LRYOKE 

c 

5 


BST=0 

c 

6 


BSY = 0 

c 

7 


BRT = 0 

c 

8 


BRY = 0 

c 

9 


ATST=0 

c 

10 


ATSY=0 

c 

11 


ATRT=0 

c 

12 


ATRY=0 

c 

13 


ATTOT=0 

c 

14 


KS AT=in 

c 

15 

C 


c 

16 

c 

STATOR TOOTH 

c 

17 

c 


c 

18 


BST=FTOTAL/ASTOTH 

c 

19 


NA = 1 

c 

20 


K=1 

c 

21 


X = BST 

c 

22 


GO TO 90 

c 

23 

10 

ATST=AT*DSS 

c 

24 

c 


c 

25 


STATOR YOKE 

c 

26 

c 


c 

27 
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20 

BSY=FPOLE/ (2. * AS YOKE) 

C 28 


NA=1 

r 29 


K=2 

C 30 


X=BSY 

C 31 


GO TO 90 

C 32 

30 

ATSY-AT+LSYOKE 

C 33 

C 


C 34 

c 

ROTOR TOOTH 

C 35 

c 


C 36 

40 

BBT=FTOTAL/ARTOTH 

C 37 


NA=31 

C 38 


K~3 

C 39 


X=BRT 

C 40 


GO TO 90 

C 41 

50 

atrt=at^dsr 

C 42 

C 


C 43 

r 

ROTOR YOKE 

C 44 

c 


C 45 

60 

BRY^FPOLE/(2. ♦APYOKE) 

C 46 


NA=31 

C 47 


K = 4 

r 48 


X=BRY 

C 49 


GO TO 90 

C 50 

70 

ATRY=AT*LRYOKE 

C 51 

C 


C 5? 

80 

attot=atag + atst+atsy+atrt+atp Y 

C 53 


RETURN 

C 54 

C 


C 55 

c 

INTERPOLATION PROCEDTIRE FOP MATERIAL CURVES 

C 56 

c 


C 57 

90 

IF (AI (NA) . LT. X) GO '^0 130 

r *^8 


NA=NA+3 

C 59 

100 

IF (AI(NA)-X) 110,120,120 

C 60 

110 

NA=NA+2 

C 61 


GO TO 100 

C 62 

120 

XX- {AI (NA) - AI (NA-2) ) / (A LOG (AX (N A + 1 ) / ( AI ( N A- 1) + 0 . 9 00 1) ) ) 

C 6 3 


Y = AI (NA) - XX^ALOG (AI (NA+ 1) ) 

C 6U 


AT = BXP ( (X- Y)/XX) 

C 65 


GO TO (10,30,50,70) ,K 

C 66 

130 

KSAT=0 

C 67 


GO TO (20,40,60,80) ,K 

C 68 


END 

C 69- 


SUBFOUTTNE SLOTS { SLTLOC , XSTYP 2 , WS X , FS X 1 , WS X2 , WSX 3 , WS X4 , WS XS , PS X , D P 1 

1 1 X,D2X,n3X,D4X,D5X,XTWD?H,CAPF^, S AREA, K, D6X, WSX6, OIA, KX,AXX, yTtvMAS p 2 


2,PHIX) p 3 

E pa 

C FOP STATOR SLOTS SLTLOC^I.O * FOR ROTOR SLOTS SLTLOC=-1.0 D S 

C P 6 

REAL N,KX D 7 

C D 8 

INTEGER XSTYPE D 9 

r D 10 

D (WA, CAREA) = { (-WA+S0KT { W A ♦ PA + 4 . S A R E A A N PH I ) ) ) / { 2. *TANPHI) D 11 

WB (P,WA) =WA+2. ♦D^TANPHI D 12 

A (W) =0.2S + W*W* ( ( (1.570B+PniX)/(COSPHI*COSPFT) ) 4-TANPHI) D 13 

C D 14 

IF (CARKA+DSX+DIX. LT. 1.0E-1F) GO TO 310 D 15 

Al=0.2S*KX+( 1.0/12.0) D 15 

A2=0.25*(KX-0. 66667) p 17 

C DIB 

GOTO ( 10, 20, 30, 90, 210, 90) , XSTYPE D 19 

C D 20 

10 WSX1=0. n 21 

D3X=0. D 22 

D4X=0. D 23 

AXX=0. D 24 

GO TO 40 n 25 

C D 26 
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20 WSXA=HSX D 27 

HSX1=0. n 2fl 

tfSX2=0- B 

D3X=0, D 30 

D4X=0. D 31 

A X X— 0 • D 32 

GO TO 100 T) 33. 

C D 34 

30 AXX = KX* (D4X/HSX 1+ (D 3X/( WSX-WSX1 ) ) ♦ (AL03 ( WSX/W.SX1) ) ) D 35 

40 WSX2=0. D 36 

WSX3=0. D 37 

WSX4=WSX D 3R 

WSX5=WSX D 39 

XTHDTH=0. D 40 

IF (DSX.GT. 1.0E-15) GO TO 50 D 41 

IF (D1X.lt. 1.0E-15) GO TO 60 D 42 

DSX:=D1X+D4X+D3X+D2X+D6X D 43 

GO TO aO D 44 

50 IF (D1X.LT.1.0E-15) GO TO 70 D 45 

GO TO aO B 46 

60 DSX=CAREA/ (WSX-2. +WSX6) +D5X+D6X+ D2X +L)3X <-D4X D 47 

70 DIX^DSX- (D4X+D3X+D2X+D6X) D 48 

80 SAREA= WSX* (D5X-D4X-D3X) 1*0. 5*(WSX1 +WSX) *D3 X+ WGX1 *04 X B 49 

IF (CAREA.T.T.I.OE-15) C A R F A= ( W3 X - 2 , * WS X 6) * ( D 1 X- D5 X) D 50 

AXX=AXX+ (D1 X*A1-D5X*A2+KX*D2X) /HSX D 51 

XTWMAG= ( (DIA + 0.6667*SLTLOC*DSX) ♦ (3. 1416/N) ) -WSX D 52 

GO TO 300 D 53 

C P 54 

90 WSXA=WSX2 B 55 

WSX = 0. D 

100 IF (KSXA ,GT. 1 -OE-15) GO TO 110 D 57 

IF (XTWDTH.lt. 1.0E-15) GO TO 110 D 58 

WSX A= ( ( 3. 1 4 16* ( DIA + 2. <'SLTLOC* (DU X^D 1 X) ) ) /N) XTRDTH D 59 

110 PHIX= (3.14 15927/N)*SLTLOC p 6Q 

GO TO 140 D 61 

12 0 XTWDTH^ ( (1 . 1 4 1 6 * (DT A + 2. *S LTLOC* ( D4X + DJX ) ) ) /N) -WGX A D 6? 

GO TO 110 B 63 

130 IF (ABS (PHIX) . LT. 1 . 0 F- 15) GOTO 120 D 6 4 

PHIX = (A BS (PMIX*0. 017453) ) *SLTL02 p 65 

XTWDTH=0. B 66 

140 IF (XSTXT^E. LT. 1 ) GO TO 150 B 67 

WSX2=WSXA B 6q 

AXX = KX* (D4 X/WSX1 ♦ (pl X/ (VSX A-WBX 1 ) ) * ( ALOG (WSX A/WGX 1) ) ) D 69 

GO TO 16P B 70 

150 WSX^WSXA n 71 

160 TANPHI = TAN (PHTX) B 72 

COSPHI=COS (PHI X) B 71 

SINPHI = SIN (PHTX) B 74 

WSX4 = WB (D2X, DSXA) D 75 

IF(DSY,gT.1.0F-15)GO'^O170 P76 

IF ( D1X ■ GT. 1. OF- 1 5) GO TO 190 D 77 

Y1 =D (HSX4-2 . *W5X6 ,CARFA/2 . ) B 79 

W = WB (D2X^-Y 1 +D5Y, WSXA) D 7D 

IF ( XSTYPE. FQ. 6) GO TO 260 B 80 

Y2 = D (W-2 .*WSX6 ,rARKA/?. ) B 81 

DSX=Y 1 +Y2 + P3X + D2X f n4X*P6X*r)5X D 8 2 

GO TO 100 B 81 

170 IF (BIX. GT. 1. OF,- 15) GO TO 200 D 84 

180 P 1 X=DSX-P4 X-Dl X-D2 X-P6 X B 85 

GO TO 200 D 86 

190 DSX=D 1 X+D4 X+D3 X+ D2 Xf D6 X n P.7 

200 IF (X5TYPE. EO. 6) GO TO 280 B 88 

WSX5=WR (02 X+D1 X, WSXA) P 89 

WSX3=WB(DSX-D3X-D4X,W.SXA) P 90 

SAREA = 0.5* (W.SXl* W8X A) * (P? Y-L4X-n IX) + 9. 5* (W.5X A^-WSX 1) *0 IX + WS/ 1*D4X P 9 1 

AXX=AXX+ ( ( 2.*KX*D2X) /( W.SXA+WSX4) ) + (0 1 X/WSX4) *A1 - { (2 .00 *P5X) / (WGX4+ D 92 

1WSX5))*A2 n 93 

XTW.1 AG= ( (DIA+2. *8LTL0C*DSX) * ( 3 . 1 4 1 6 / N ) ) - VS Y 3 D 04 

IF (CAPEA.GT. 1.0E- 15) GO TO 300 D 95 

A3=ni X* (WSX5+ WSX4-4 . * WSX6 ) D 96 

CAR EA = C. 5*A3-P5X *SpPT ( ( '■;SX5-2.* V SX6) **2-A3* NP!! r - { D5 X*'^’ A NPF D **2) 0 0 7 

GO TO 300 p 9H 

C P 99 

210 WSX=0 p 100 
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u tJ u 


HSX2=0. D 101 

WSX4=0. D 102 

XTWDTH = 0. D 103 

D2X~0. D 104 

D3X=0. D 105 

D5X=0, D 106 

WSX5=0. D 107 

WSX6=0, D 108 

IF (DSX,GT.1.0E-15) GO TO 240 D 109 

IF (D1X.lt. 1,0E-15) GO TO 220 D 110 

DSX=D1X+2. ♦D6X+D4X D 111 

GO TO 250 D 112 

220 DSX= (S0RT(4.*CAREA/3. 1416) ) +D4X+2.^D6X D 113 

230 D1X=DSX-2.*D6X-D4X D 114 

GO TO 250 D 115 

240 IF (D1X.LT.1.0E-15) GO TO 230 D 116 

IF (ABS (DlX + D4X+2.*D6X-DSX),Gr. 3,001) GO TO 310 D 117 

250 IF (WSX2.LT. 1.0E- 15) WSX2=DSX-D4X D 11R 

SAREA=0.7854*WSX2*WSX2+WSX1*D4X D 119 

IF (CAPEA.LT. 1.0E-15) C AH EA ^0 . 78 54* ( WSX2- 2 . *06 X) D 120 

AXX=(0.625+D4X/WSX1)*KX D 121 

XTWMAG= (DIA* (2, * D4X + 1 . 3 3 3* WSX2) *SLTLOC) * (3. 1416/N) -0. 94*WSX2 D 122 

GO TO 340 D 123 

C D 124 

260 W1=W-2,*D6X D 125 

AS=A (HI ) D 126 

IF (2, *AS/CAPEA. GT. 1 .00) GO TO 310 D 127 

270 AE=0.5*CAREA-AS D 128 

Y2=D(W1,AR) D 129 

W2=WB(Y2,W1) D 130 

AS = A(H2) D 131 

IF (ABS (2. * (AR* AS) /CAREA-1.0) . GE. 0. 00 1) GO TO 270 D 132 

DSX=D4X + D3X + D2X+D6X+I H- Y2 +D5X+ ( W 2 /2 . ) * ( 1 . /C OSPH I + A N PH T) D 133 

D1X=DSX-D4X-D3X-D2X-D6X D 134 

280 WSX3= (2. * (DSX-D4X-D3X) *SINPHI*WSX2*COSPHI) / (SINPRI* 1 . 0) D 135 

WSX5=HSX3 D 136 

SAREA = A (WSX3*COSPHI) <-0. 5* (WSX 3* C OS PHI + WSX 2) * (DSX-0. 5*WSX3* { 1 . 0+STN P 137 

1PHI) -D4X-D3X) +0.5*D3X+(WSX2+WSX1) +WSX1+D4X D 138 

AXX=AXX+ ( (2. *KX*D2X) / (WSX A+W5X4) ) + (01X/WSX4) *A 1“ ( (2.00+D5X) / (WSX4+ D 139 

1WSX5) ) *A2 D 140 

XTHMAG= ( (DIA+2. *SLTLOC* (D4X + D3X) ) ♦ (3. 14 16/N) ) -WSX2 D 141 

IF (CAREA.GT.1.0E-15) GO TO 300 D 142 

Y1=0.5*DSX D 143 

W = 0.5* (WSX4+WSX5) D 144 

CAREA2-=1000. D 145 

290 CAREA=A (WSX3*C0SPHI-2. *D6X) +0. 5* (WSX3*COSPKI + H“4. *D6X) * (OSX-0.5* WS D 146 

1X3* ( 1. 0+STN PH I) - D4X-D3X-D2X-Y1-D5X) + ( 0 . 5* { Y 1 + D5 X) * ( W+ NS X 4- 4 . *WSX 6) D 1 47 

2)- (W-2. *WSX6) *D5X D 148 

IF (ABS (CAREA2/CARSA-1 . 0) . LT.O. 00 1) GO TO 300 D 149 

Yl=D(WSX4“2.*HSX6,CAHEA/2.) D 150 

W=WB (Y 1 +D5X+D2X, WSX2) D 151 

CAREA2=CAREA n 152 

GO TO 290 D 153 

C r 154 

300 IF ( ABS ( D1X+D2X + D3X+D4X + D6X-D SX) .LT. 0. 001 ) GO TO 340 D 155 

C n 1 56 

310 IF (SLTLOC.LT.O.) WRITE (6,320) ' D 157 

IF (SLTLOC.GT.O. ) WRITE (6,330) D 158 

320 FORMAT ( 1 HK , 4 1 HINSUFFIC TENT OR INCORRECT ROTOR SLO-^ DATA) D 159 

330 FORMAT ( 1 H K,4 2 HI NS DFFIC I E NT OR INCORRECT STATOR SLOT DATA) D 160 

340 RETURN D 161 

END D 162- 


SOBROUTINE WDGFCT (PB A , P , OS , DF , PC, PF, WDGPCH) E 1 

E 2 

PITCH FACTOR CALCULATION E 3 

E 4 

YY=FLOAT (IFIX ( { (QS/P) *WDGPCH) +0.01)) E 5 

IF (ABS (TY-QS/P*WDGPCH) .GT. 1 .OE-2) WRITE (6,10) WDGPCH E 6 


10 FORMAT ( 1HK,F5. 3,22H PITCH IS NOT POSSIBLE) 

PF=SIN (YY*1 ,571/ (QS/P) ) 

C 

C DISTRIBUTION FACTOR CALCULATIONS 

C 

IPX=IFIX(P+0. 1) 

IQQ=IFIX {QS + 0. 1) 

IC=IFIX(PC+0.1) 

IPN=3 

PN=3. 

QN=QS/(3,*P) 

C 

C CHECK IF WINDING HAS INTEGRAL NO. OF SLOTS PER POLE PER PHASE 

C 

D = 1.0 

IF (PBA. GT. 61. 0) D^2.0 

IZY=IPX*IPN 

IDM = 0 

20 IDM=IDM+IZY 

IF (IQO-IDM) 40,30,20 
C 

C CALCULATE DISTRIBUTION FACTOR FOR INTEGRAL SLOT WINDING 

C 

30 DF=SIN ( 1. 571*D/PN) /(QN*D*STN ( 1 . 5 7 1 / (PN* QN ) ) ) 

GO TO 90 
C 

C REDUCE THE FRACTION IQQ/IZY TO LOWEST TERMS 

C 

40 II0Q=IQQ 

1 = 2 

50 IF ( (IZY/I) *I . EQ.I2Y. AND. (TIOQ/X) *T. Eg. IIOQ) GO TO 60 
IF (I. GT. IZY) GO TO 70 
J=I* 1 
GO TO 50 
60 IZY=IZY/I 

IIQg=IIOQ/T 
GO TO 50 
C 

C CALCULATE DISTRIBUTION FACTOR FQH FRACTIONAL SLOT WINDING 

C 

70 FNQ=IIQQ 

DF=STN (1 .571*D/PN) / (FNQ*D^SIN ( 1 . 57 1 / ( FNg* PN) ) ) 

IF ( (IZY/3) ♦ 3. EO. T2Y) WRITE (6,80) 

IF ( (IPX/IZY) *IZY.NE.IPX) WRITE (6,80) 

80 FORMAT ( IHK , 40H1 M PRO PER F R ACTION A L- SLOT WINDING IS USED) 

C 

C CHECK IF SPECIFIED NUMBER OF P\RALLEL CIRCUITS ARE POSSIBLE 

C 

IPX= IPX/IZY 

90 IF ( (IPX/IC) *IC. SQ. I PX) GOTO 110 
WRITE (6,100) IC 

100 FORMAT (1HK,I2,3SH PARALLEL CIRCUITS ARE NOT POSSIBLE) 

110 RETURN 
END 


E 7 
E 8 
E 9 
E 10 
E 11 
E 12 
E 13 
F 14 
E 15 
E 16 
E 17 
E 18 
E 19 
E 20 
E 21 
E 22 
E 23 
E 24 
E 25 
E 26 
E 27 
E 28 
E 29 
E 30 
E 31 
E 32 
E 33 
E 34 
E 35 
E 36 
E 37 
E 38 
E 39 
P. 40 
E 41 
F 42 
E 43 
P 44 
E 45 
E 46 
E 47 
F 48 
E 49 
E 50 
E 51 
E 52 
E 53 
E 54 
E 55 
F 56 
E 57 
E 58 
F 59 
E 60- 


SUBROUTINE CMBNTN (QS,NB,P) 
C 

REAL NB 

DIMENSION L (100) 

C 

X=1.0E-15 

K=0 

F=NB 

10 D=ABS(QS-F) 

?i = 1 

DO 20 1=1, 1000 
A=3,*FLOAT {I)*P 
IF (ABS (D-A) .LT.X) GO TO 40 
IF (A.GT.D) go to 30 


F 1 
F 2 
F 3 
F 4 
F 5 
F 6 
F 7 
F 8 
F 9 
F 10 
F 11 
P 12 
F 13 
F 14 
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20 CONTINUE F 15 

30 IF (ABS(D-P) .LT. X) GO TO 40 p 16 

IF {ABS (P-FLOAT (IFIX (F/P+0.0001) ) ) . LT.X) GO TO 40 F 17 

M=2 F 18 

IF (F.GT.QS+P/2. ) GO TO 40 F 19 

M=3 F 20 

IF (ABS {D-P/2.) -LT. X) GO TO 40 F 21 

IF (ABS (QS-F) .LT.X) GO TO 40 F 22 

IF (ABS (D-1 .) . LT. X) GO TO 40 F 23 

IF (ABS (D-2.) -LT. X) GO TO 40 F 24 

IF (ABS {D“P+1 - ) . LT. X) GO TO 40 F 25 

IF (ABS (D-P-1.) .LT.X) GO TO 40 F 26 

IF (ABS (D-P-2. ) . LT. X) GO TO 40 F 27 

IF (ABS (D-P+2.) .LT.X) GO TO 40 F 28 

IF (K.EQ.O) GO TO 150 F 29 

11=11+1 F 30 

L (II ) =IFIX (F + 0.0 1 ) F 31 

GO TO 110 F 32 

40 IF (K.GT.O) GO TO 110 F 

K=1 F 34 

F = FLOAT (IFIX (0.60^0S) ) F 35 

11=0 F 36 

C F 37 

GO TO (50,70,90) F 38 

C F 39 

50 WRITE (6,60) F 40 

60 FORMAT ( 1 HK , 82 HROTO R-5TA TOP SLOT COMBINATION MAY PRODUCE UNDESTRAB F 41 

1LE TOSOUE-SPBED CHARACTERISTICS) F U2 

GO TO 110 F 43 

C F 44 

70 WRITE (6,60) F 45 

FF=(P/2.)/(P/2.+QS) F 46 

WRITE (6,80) FF F 47 

80 FORMAT (1H ,19HMINIMI2E BY S KEWI NG, F6. 3, 30 H TIMES ROTOR CIRCUMFERE F 48 

1NCE, OR) F 49 

GO TO 1 10 F 50 

C F 51 

90 WRITE (6,100) F 52 

100 FORMAT (1HK,61HP0T0R-STAT0P SLOT COMBINATION MAY PRODtlCF NOISE AND F 53 

1 VIBRATION) F 54 

C F 55 

110 F=F+1.0 F 56 

IF (F. LE. 1 .4*QS) GO TO 10 F 57 

IF (II. LT. 1) GO TO 1 30 F 58 

C F 59 

WRITE (6,120) (L (I) ,1=1 ,11) F 60 

120 FORMAT (1H ,52HCHANGE NUMBER OF ROTOR SLOTS TO ONE OF THE FOLLOHTN F 61 

1G/(10I6)) F 62 

GO TO 150 F 63 

C F 64 

130 WRITE (6,140) F 65 

140 FORMAT (1H ,29HCHANGE NUMBER OF STATOR SLOTS) F 66 

C F 67 

150 RETURN F 68 

END F 69- 


50 



APPENDIX D 


ERROR MESSAGES 

This appendix lists the various error messages that may result diming program ex- 
ecution. For each error message the subroutine from which the message originated is 
identified and the probable cause of the error is suggested. The purpose of these error 
messages is only to warn and to inform. In no case is program execution terminated. 
This information is provided in the following table: 


Num- 

ber 

Error message 

Responsible 

subroutine 

Explanation 

1 

CORE LOSS DATA IS NOT GIVEN 
AT SPECIFIED STATOR LAMINATION 
THICKNESS 

USE DATA FOR xx. xxx LAMINATIONS 

INDMTR 

(1) The stator material deck does not 
contain core- loss data card 
($FELOSS) for lamination thickness 
within 0. 0005 in. of lamination thick- 
ness specified on the motor design 
deck data card $STATOR. The pro- 
gram will use the best available 
core -loss data. 

(2) Core-loss data may have been 
omitted entirely. 

2 

INSUFFICIENT STATOR SLOT DATA, 
SPACE FACTOR OF 0. 70 ASSUMED 

INDMTR 

DSS, DIS, SCAREA, and CSRATO 
are all less than 1. OE-15. The pro- 
gram assumes a value of CSRATO = 
0.70. 

3 

SHUNT RESISTANCE RO FAILED TO CONVERGE 

INDMTR 

The iteration for RO in the no-load 
magnetic calculations did not converge 
after 10 iterations. This generally 
means that the magnetic flux path is 
saturated or nearly saturated. 

4 

MAGNETIZING CURRENT FAILED TO CONVERGE 

INDMTR 

The iteration for magnetizing current 
and XO in the no-load magnetic cal- 
culations did not converge after 15 
iterations. This generally implies 
that the motor is magnetically satu- 
rated or nearly saturated or that the 
material has square -loop character- 
istics with the flux density falling 
near the knee of the curve. 
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Num- 

ber 

Error message 

Responsible 

subroutine 

Explanation 

5 

MACHINE SATURATED 

INDMTR 

One or more parts of the magnetic 
circuit of the motor saturated at no 
load. In order to determine which 
part or parts, compare the computed 
flux densities with the maximum flux 
density for the appropriate material. 
The ampere-turn drop across any 
part of the magnetic circuit that sat- 
urated is assumed to be zero. 

6 

INSUFFICIENT DATA TO SCALE WINDAGE LOSS 

INDMTR 

(1) One or more of the following var- 
iables is very small or zero: DIAREF, 
LREF, RPMREF, GAPREF. All of 
these variables must be defined to 
permit scaling of windage loss. 

(2) The variable WL is very small 
or zero. The synchronous windage 
loss will be assumed to be zero. 

7 

F + W TORQUE EXCEEDS AVAILABLE SHAFT 
TORQUE AT xxx.xx PERCENT SLIP 

INDMTR 

1 

This message is printed if the total 
electromagnetic shaft torque computed 
in subroutine CIRCT is less than the 
! computed windage torque at the speci- 
fied value of slip. If this error oc- 
curs for values of slip greater than 
15 percent, equivalent circuit analy- 
sis is terminated and the program 
proceeds to plot the torque-speed 
curve. For values of slip below 
15 percent the program continues to 
increment slip in the normal manner. 

8 

INSUFFICIENT OR INCORRECT ROTOR SLOT 
DATA 

SLOTS 

(1) SB, DSR, and DIR are all less 
than 1. OE-15. At least one of these 
variables must be read in. 

(2) For slot type 6 only: if area AR 
(fig. 14) becomes negative, this mes- 
sage is printed. In order to eliminate 
this problem, make the slot narrower 
and deeper. 
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Num- 

ber 

Error message 

Responsible 

subroutine 

Explanation 

9 

INSUFFICIENT OR INCORRECT STATOR SLOT 
DATA 

SLOTS 

(1) DSS, DIS, and SCAREA are all 
less than 1. OE-15 at the time subrou- 
tine SLOTS is called to compute stator 
slot dimensions. In general, this im- 
plies one of the following: 

CSS is zero or negative 
ASTRND is zero or negative 
AWG is not between 1 and 40 inclu- 
sive 

STRNDS is zero or negative. 

(2) For slot type 6 only: see error 
message 8. 

10 

.XXX PITCH IS NOT POSSIBLE 

WDGFCT 

This message is printed if the number 
of stator slots per pole times the 
stator winding pitch WDGPCH is not 
within 0. 01 of an integer value. 

11 

XX PARALLEL CIRCUITS ARE NOT POSSIBLE 

WDGFCT 

This message is printed for either 
fractional or integral slot windings. 

It means that a balanced, three-phase 
winding is not possible with the num- 
ber of parallel circuits specified in 
the input data. 

12 

IMPROPER FRACTIONAL SLOT WINDING 
IS USED 

WDGFCT 

(1) The denominator of the slots per 
pole (reduced to lowest common de- 
nominator) is not divisible by 3. 

(2) The number of poles is not divis- 
ible by the denominator of the slots 
per pole per phase (reduced to lowest 
terms). 

13 

ROTOR- STATOR SLOT COMBINATION MAY 
PRODUCE UNDESIRABLE TORQUE-SPEED 
CHARACTERISTICS 

CMBNTN 

See reference 3 (pp. 317-320) and 
reference 5. 

14 

ROTOR-STATOR SLOT COMBINATION MAY 
PRODUCE NOISE AND VIBRATION 

CMBNTN 

See reference 3 (pp. 317-320) and 
reference 5. 

15 

MINIMIZE BY SKEWING x. xxx TIMES 
ROTOR CIRCUMFERENCE, OR 

CMBNTN 

This message can only follow error 
message 13 and is always followed by 
message 16 or 17. It states that the 
undesirable effects referred to in 
error message 13 can be reduced or 
eliminated by skewing the specified 
amount. 
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Num- 

ber 

Krror message 

Responsible 

subroutine 

Explanation 

16 

CHANGE NUMBER OF ROTOR SLOTS TO 
ONE OF THE FOLLOWING 

XXX XXX XXX 

CMBNTN 

This error message follows error 
message 13, 14, or 15 but never in 
conjunction with message 17, It lists 
the number of rotor slots that may be 
used without incurring the problems 




referenced in message 13 or 14. 

17 

CHANGE NUMBER OF STATOR SLOTS 

CMBNTN 

This message follows error message 
13, 14, or 15 but never in conjunction 
with message 16. It is displayed only 
if it is not possible to find a number 
of rotor slots that will eliminate the 
problems referenced in error mes- 
sage 13 or 14. 
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APPENDIX E 


ALPHABETIC FORTRAN SYMBOL LIST 

An alphabetic FORTRAN symbol list is given for the main program and each sub- 
routine. The S 3 mibol list for the main program is given first. This list is complete, 
showing every symbol used in the main program. The symbol lists for the subroutines 
follow in this order: SLOTS, CMBNTN, WDGFCT, MAGNET, and CIRCT. The symbol 
lists for the subroutines list only those FORTRAN variables that do not appear in the 
main program or those which, if they do appear in the main program, have a definition 
different from that in the main program . 

Where S 3 unbols define stator or rotor slot dimensions, further clarification may be 
obtained by referring to figures 4, 12, and 14. Figure 4 shows all slot dimensions that 
are needed to calculate the slot permeance ratio. Figure 12 shows all slot dimensions 
that are allowable input. Figure 14 shows those slot dimensions that are not shown in 
either of the other two figures. 


AI 

AIRGAP 

ARTOTH 

ARYOKE 

ASTOTH 

ASTRND 

ASYOKE 

ATAG 

ATRT 

ATRY 

ATST 

ATSY 

ATTOT 

AWG 


Main Program 

coordinates of points on rotor and stator material magnetization curves 
NAMELIST name 

2 

cross-sectional area of rotor teeth (used in magnetic calculations), in. 

2 

cross-sectional area of rotor yoke (used in magnetic calculations), in. 

2 

cross-sectional area of stator teeth (used in magnetic calculations), in. 

2 

cross-sectional area of stator strand, in. 

2 

cross-sectional area of stator yoke (used in magnetic calculations), in. 

ampere-turns across airgap, ampere-turns 

ampere-turns across rotor tooth, ampere-turns 

ampere-turns across rotor yoke, ampere-turns 

ampere-turns across stator tooth, ampere-turns 

ampere-turns across stator yoke, ampere-turns 

total ampere-turn drop, ampere-turns 

strand size of stator winding (American Wire Gage) 
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AXR 

AXS 

AY 

AY 

B 

BG 

BK 

BLANK 

BR 

BRT 

BRY 

BST 

BSY 

CO 

Cl 

C2 

C3 

C4 

CALPHA 

CCR 

CCS 

cm 

CIRCT 

CLOSS 

CMBNTN 

CSRATO 

CSS 


rotor slot leakage permeance ratio 
stator slot leakage permeance ratio 
length of one end-turn, in. 

multiplier in slot and end- connection reactance calculations and in rotor 
resistance calculations 

armature coil extension, in. 

2 

average airgap flux density, kilolines/in. 

2 

flux density at which WFE and WCORE are specified, kilolines/in. 
storage location for storing a BCD blank 

spacing between end-ring and rotor laminations (ref. 3, p. 336, fig. 199), 
in. 

2 

flux density in rotor tooth, kilolines/in. 

2 

flux density in rotor yoke, kilolines/in. 

2 

flux density in stator tooth, kilolines/in. 

2 

flux density in stator yoke, kilolines/in. 
coefficient of viscosity polynomial (see VSCSTY) 
coefficient of viscosity polynomial (see VSCSTY) 
coefficient of viscosity polynomial (see VSCSTY) 
coefficient of viscosity polynomial (see VSCSTY) 
coefficient of viscosity polynomial (see VSCSTY) 
cosine (alpha) (ref. 3, p. 209, fig. 135) 

Carter coefficient (rotor) 

Carter coefficient (stator) 
common block name 
subroutine name 
array containing core- loss data 
subroutine name 

space factor (=CSS* SS/SCAREA) 
number of conductors per stator slot 
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CURDEN 

D 

DIR 

DIS 

D2R 

D2S 

D3R 

D3S 

D4R 

D4S 

D5S 

D6R 

D6S 

DBRS 

DBS 

DC 

DEKTYP 


DELTAS 

DERI 

DER2 

DIAREF 

DIFF 

DIFFl 

DIR 

DNSTY 


2 

current density in armature, A/in. 

stator bore diameter, in. 

overall conductor depth in rotor slot, in. 

overall conductor depth in stator slot, in. 

rotor slot dimension, in. 

stator slot dimension, in. 

rotor slot dimension, in. 

stator slot dimension, in. 

rotor slot dimension (slot-opening depth), in. 

stator slot dimension (slot- opening depth), in. 

stator slot dimension, in. 

rotor slot dimension, in. 

stator slot dimension, in. 

depth below rotor slot, in. 

depth below stator slot, in. 

distance between center of end-ring and center of stator slot (ref. 3, p. 336, 
fig. 199), in. 

character in card column 1 of first card following each motor design deck 
if DEKTYP := M, it marks start of new data set; if DEKTYP = BLANK it 
marks start of new motor design data deck 

increment by which S is increased, percent 
end- ring outside diameter, in. 
end-ring inside diameter, in. 
reference diameter for scaling windage loss, in. 
smallest of all values of DIFFl calculated, in. 

difference between stator lamination thickness and lamination thickness 
specified on $FELOSS data card, in. 

rotor lamination inside diameter, in. 

array containing density values for various rotor and stator winding material 

Q 

possibilities, Ib/in. 
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DOS 

DR 

DSR 

DSS 

EFF 

ENDTRN 

F 

FI 

FCORE 

FELOSS 

FLDNME 

FPOLE 

FTOTAL 

FW 

FWl 

G 

GAPREF 

GE 

HP 

I 

II 
lA 

IBAR 

ICNTl 

ICNT2 

IMAG 

IMAG2 


stator lamination outside diameter, in. 

rotor lamination outside diameter, in. 

rotor slot depth, in. 

stator slot depth, in. 

efficiency, percent 

axial length of end turn, in. 

frequency of line voltage, Hz 

part of horizontal extension of armature winding (ref. 3, p. 209, fig. 135), 
in. 

frequency at which WCORE is given, Hz 
NAMELIST name 

name of fluid in motor cavity (must be limited to six characters or less) 

flux per pole, kilolines 

total flux, kiloUnes 

windage loss at rotor speed (rpm), W 

windage loss at synchronous speed, W 

airgap, in. 

reference gap for scaling windage loss, in. 

effective airgap, in. 

shaft power, hp 

subscript or index 

line current, A 

subscript or index 

rms current in one rotor bar, A 

counts number of iterations on RO during no-load magnetic calculations 

counts number of iterations on magnetizing current during no-load magnetic 
calculations 

magnetizing current, A 
magnetizing current, A 
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INITL 

J 

JBAR 

JRING 

KDS 

KOBE 

KPS 

KRING 

KS 

KSAT 

KT 

L 

LARM 

LAST 

LB 

LREF 

LRYOKE 

LS 

LSYOKE 

LT 

LTOTAL 

LTR 

LTS 

MAG 

MAGNET 

MATDEK 


common block name 
subscript or index 

2 

current density in rotor bar, A/in. 

2 

current density in end ring, A/in. 
distribution factor for stator winding 
input to plotting routine PLOTXY 
pitch factor for stator winding 

correction factor for end-ring resistance (ref. 3, p. 334, fig. 194; and 
ref. 4) 

slot leakage pitch factor (ref. 2, p. 185, fig. 7.3) 

saturation indicator 

index 

stator core length, in. 

total length of wire of armature winding, ft 

logical variable - LAST=. TRUE. - indicates last core- loss data card has 
been read 

length of rotor bar (including portion inserted in end-ring), in. 
reference length for scaling windage loss, in. 
length of flux path through rotor yoke, in. 

length of one armature conductor (half of armature coil length), in. 
length of flux path through stator yoke, in. 

thickness of laminations at which core- loss data are given in material deck, 
in. 

overall axial armature length (2. *ENDTRN + L), in. 

thickness of rotor laminations, in. 

thickness of stator laminations, in. 

common block name 

subroutine name 

alphabetic constant (defined to be character ”M” in a data statement) 
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N 

NAME 

NB 

NCARDS 

NSYNCH 

P 

PC 

PF 

PFLUID 

PHASE 

PHIR 

PHB 

PIN 

PLOTXY 

POUT 

PP 

PREF 

QS 

RO 

ROOLD 

R1 

R2 

R2BAR 

R2RING 

RATING 


number of stator conductors in series per phase (2* (number of stator turns 
in series per phase)) 

subscripted array containing information in columns 3 to 80 of first card fol- 
lowing each motor design deck 

number of rotor bars (equal to number of rotor slots) 

number of core- loss data cards ($FELOSS) read in (last card ($FELOSS 
LAST=. TRUE. $) is not counted) 

synchronous speed of motor, rpm 

number of poles 

number of parallel circuits 

power factor 

pressure of fluid in airgap, psi 

if PHASE equals BCD BLANK, PF is lagging; if PHASE equals * , PF is 
leading 

one-half of angle at which rotor slot sides diverge, deg 
one-half of angle at which stator slot sides diverge, deg 
power input to motor, W 
subroutine name 

output power available at motor shaft, W 
input to plotting routine PLOTXY 

reference pressure of fluid in airgap used for scaling windage loss, psi 
number of stator slots 

shunt resistance of equivalent circuit, ohms 

value of RO calculated during previous iteration pass, ohms 

armature resistance, ohms 

rotor resistance referred to stator winding, ohms 
component of R2 attributable to rotor bars, ohms 
component of R2 attributable to end rings, ohms 
NAMELIST name 
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RATIO 

RESETl 

RESET2 

RMAT 

ROTOR 

RPM 

RPMREF 

RRSTVY 

RSAREA 

RSLOTS 

RSTVTY 

RSTYPE 

RTRWDG 

RTWDTH 

RTWMAG 

RWMAT 

S 

SALPHA 

SB 

SCAREA 


SER 

SFR 

SFS 

SKEW 


WSS4/WSS3 for trapezoidal stator slot; WSR4/WSR3 for trapezoidal rotor 
slot; DER2/DER1 for rotor-winding end-ring 

array made equivalent to common block INTTL 

array made equivalent to first seven entries in common block CIR 

array containing description of rotor lamination material 

NAMELIST name 

rotor speed at slip S, rpm 

reference RPM for scaling windage loss, rpm 

resistivity of rotor winding material at temperature TRW, pin. -ohm 
2 

rotor slot area, in. 

NAMELIST name 

array containing resistivity values for various rotor and stator winding 
materials at 20° C, pin. -ohm 

rotor slot type 

NAMELIST name 

rotor tooth width (if constant), in. 

rotor tooth width used in magnetic calculations, in. 

code for rotor winding material: 1 for aluminum; 2 for brass; 3 for copper 

clearance between armature coils at end turns (ref. 2, p. 309, table 26; and 
p. 209, fig. 135), in. 

sin (ALPHA) (ref. 2, p. 209, fig. 135) 

2 

cross-sectional area of rotor bar, in. 

slot area remaining after subtracting, from total slot area, slot opening and 

approximate areas occupied by slot liners, separators, wedges, etc. 

2 

(shaded area in fig. 12), in. 

2 

end-ring cross-sectional area, in. 
rotor lamination stacking factor 
stator lamination stacking factor 

skew of rotor slots measured along rotor circumference, in. 
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SLIP 

SLOPE 

SLOTS 

SMAT 

SMAX 

SOLD 


SPITCH 

ss 

SSAREA 

SSLOTS 

SSTYPE 

STATOR 

STRNDS 

STRWDG 

STWDTH 

STWMAG 

SWMAT 

T 

TIR 
TIS 
TER 
T FLUID 
TITLE 

TMPCF 


array containing values of slip at which motor performance is calculated, 
percent 

slope of core-loss- against- frequency curve (for constant flux density) on 
log- log graphs, measured at frequency FCORE and flux density BK 

subroutine name 

array containing description of stator lamination material 

maximum value of S for which motor performance is calculated, percent 

previous value of S at which motor performance was calculated (used to 
calculate S at rated torque and to resume calculations at proper value 
of S following calculations at rated torque), percent 

stator winding pitch expressed as a decimal fraction, per unit 

2 

cross-sectional area of stator conductor, in. 

2 

total area of stator slot, in. 

NAMELIST name 
stator slot type 
NAMELIST name 

number of strands per armature conductor 
NAMELIST name 

stator tooth width (if constant), in. 

stator tooth width used in magnetic calculations, in. 

code for stator winding material: 1 for aluminum; 2 for brass; 3 for copper 

shaft torque at slip S, in. - lb 

rotor slot pitch at airgap, in. 

stator slot pitch at airgap, in. 

end-ring thickness, in. 

temperature of fluid in motor cavity, °C 

array which contains name or description of design to be analyzed, used to 
print heading on output listing 

array containii^ temperature coefficients of resistivity for various possible 
rotor and stator winding materials, per °C 
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TOLD 

TORQUE 

TRATED 

TREE 

TRW 

TSW 

VI 

V2 

VSCFLD 

VSCREF 

VSCSTY 


WO 

W1 

W2 

WAREA 

WARM 

WEAR 

WCORE 

WDGFCT 

WEIGHT 

WFE 

WL 

WNDAGE 

WRING 

WROT 


value of T at previous value of S, in. -lb 

array containing values of T corresponding to values of S stored in array 
SLIP, in. -lb 

rated torque, in. - lb 

reference temperature for scaling windage loss, °C 

temperature of rotor winding, °C 

temperature of stator winding, °C 

line-to-neutral voltage, rms volts 

airgap voltage, rms volts 

viscosity of fluid in motor cavity, Ibm/ft-sec 

reference viscosity for scaling windage loss, Ibm/ft-sec 

arithmetic statement function, VSCSTY = CO + Cl* T + C2* T* * 2 + C3* T* * 3 
+ C4*T**4, where VSCSTY is fluid viscosity in Ibm/ft-sec and T is 
fluid temperature in °C; CO to C4 are program input 

core loss, W 

losses in armature winding, W 
losses in rotor winding, W 

2 

array containing cross-sectional areas of standard wire gages, in. 
weight of armature (exclusive of insulation), lb 
power loss in one rotor bar, W 

core loss for stator laminations at frequency FCORE and at flux density 
BK, W/lb 

subroutine name 

total electromagnetic weight, lb 

core loss for stator laminations at frequency F and at flux density BK, 
W/lb 

windage loss at reference conditions, W 
NAMELIST name 
loss per end-ring, W 
rotor iron weight, lb 
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WRWNDG 

WSR 

WSRl 

WSR2 

WSR3 

WSR 4 

WSR5 

WSR6 

WSS 

WSSl 

WSS2 

WSS3 

WSS4 

WSS5 

WSS6 

WSTAT 

WSTOTH 

WSYOKE 

XO 

XOAG 

XI 

X2 

XLGND 

XP 

XRE 

XRS 

XRZ 

XSE 


weight of rotor winding, lb 
rotor slot width (if constant), in. 

width of rotor slot opening (for partially closed slot), in. 

rotor slot dimension, in. 

rotor slot dimension, in. 

rotor slot dimension, in. 

rotor slot dimension, in. 

rotor slot dimension, in. 

stator slot width (if constant), in. 

width of stator slot opening (for partially closed slot), in. 

stator slot dimension, in. 

stator slot dimension, in. 

stator slot dimension, in. 

stator slot dimension, in. 

stator slot dimension, in. 

stator iron weight, lb 

weight of stator teeth, lb 

weight of stator yoke (back iron), lb 

magnetizing reactance, ohms 

magnetizing reactance of airgap only, ohms 

armature leakage reactance, ohms 

rotor leakage reactance referred to stator winding, ohms 

array containing legend printed to left of slip-torque plot 

peripheral airgap leakage reactance, ohms 

rotor end- turn leakage reactance, ohms 

rotor slot leakage reactance, ohms 

rotor zigzag reactance, ohms 

stator end-turn leakage reactance, ohms 
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XSK one- half of total skew reactance, ohms 

XSS stator slot leakage reactance, ohms 

XSZ stator zigzag reactance, ohms 

XX index used dxiring no-load magnetic calculations: 1. 0 if XO is to be calcu- 

lated; 0. if XO was read in 

XY index used during no-load magnetic calculations: 1. 0 if RO is to be calcu- 

lated; 0. if RO was read in 

XZ multiplier for zigzag reactances 


Subroutine CIRCT 

Definitions of those variables that are not listed are the same as in the main pro- 
gram. 

A real part of various complex variables 

B imaginary part of various complex variables 

C constant (C = 2. 5) 

D determinant of coefficients of circuit equations 

Fl complex input voltage to equivalent circuit (Une-to-neutral input voltage to 

motor), rms 

F2 complex voltage across shunt branch of equivalent circuit, rms 

12 current through Z2, A 

lA complex current through Zl, A 

IB complex current through Z2, A 

IC complex current through ZO, A 

STAR storage location storing BCD character * 

ZO impedance of shunt branch of equivalent circuit, ohms 

Zl stator impedance, ohms 

Z2 rotor impedance referred to stator, ohms 
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Subroutine MAGNET 


Definitions of those variables that are not listed are the same as in the main pro- 
gram. 

AT ampere-turn drop across various sections of magnetic circuit, ampere-turns 
NA subscript 
K index 

X flux density at which AT is found by interpolation between points on magnetization 

2 

curve, kilolines/in. 

XX slope of magnetization curve at flux density X 
Y used in interpolation procedure for AT 


Subroutine SLOTS 


Definitions of those variables that are not listed are the same as in the main pro- 
gram. 


A 

A1 

A2 

AR 

AS 

AXX 

CAREA 

CAREA2 

COSPHI 

D 

DIX 

D2X 

D3X 


arithmetic function 

constant used in slot permeance ratio calculations 
constant used in slot permeance ratio calculations 

slot area (fig. 14) needed for intermediate calculations for slot type 6 only, 

in. 2 

slot area (fig. 14) needed for intermediate calculations for slot type 6 only, 

in. 2 

slot leakage permeance ratio 

slot area remaining after subtractii^ slot opening, slot liners, separator, 

. . 2 
etc. , in. 

value of CAREA during a previous iteration pass (used with slot t3rpe 6 only) 
cos (phi) 

arithmetic function 
slot dimension, in. 
slot dimension, in. 
slot dimension, in. 
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Subroutine SLOTS 


D4X 

D5X 

D6X 

DIA 

DSX 

KX 

N 

PHIX 

SAREA 

SINPffl 

SLOTS 

SLTLOC 

TANPHI 

W 

W1 

W2 

WA 

WB 

wsx 

WSXl 

WSX2 

WSX3 

WSX4 

WSX5 

WSX6 

WSXA 


slot dimension, in. 
slot dimension, in. 
slot dimension, in. 

rotor outside diameter if SLTLOC = -1. 0; stator inside diameter if 
SLTLOC = 1.0, in. 

slot dimension, in. 

equals 1. 0 for rotor slots; equals slot leakage pitch factor for stator slots 
(ref. 2, p. 185, fig. 7.3) 

number of slots 

one-half of angle at which slot sides diverge (PHIX is negative for rotor 
slots, positive for stator slots), rad 
2 

total slot area, in. 
sin (phi) 

subroutine name 

indicates slot location: 1. 0 for stator slots; 1. 0 for rotor slots 
tan (phi) 

slot dimension, in. 
slot dimension, in. 
slot dimension, in. 

dummy variable used in arithmetic function definition 

arithmetic function 

slot dimension, in. 

slot dimension, in. 

slot dimension, in. 

slot dimension, in. 

slot dimension, in. 

slot dimension, in. 

slot dimension, in. 

equals WSX for slot t3q)e 2; equals WSX2 for slot types 4 and 6, in. 
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Subroutine SLOTS 


XSTYPE 

XTWDTH 

XTWMAG 

Y1 

Y2 


slot t 3 ?pe 

tooth width (for slot types 2, 4, and 6 only), in. 

average tooth width used in magnetic calculations in subroutine MAGNET, in. 
slot dimension, in. 
slot dimension, in. 


Subroutine WDGFCT 

Definitions of those variables that are not listed are the same as in the main pro- 
gram. 

D constant: 1. 0 for windings with phase belt less than 60°; 2. 0 for windings 

with phase belt greater than 60° 

DF distribution factor 

FNQ real variable equal to IIQQ after fraction ’’slots per pole per phase” has 

been reduced to lowest terms 

I integer that is tested to see if it is a common divisor of fraction ’’slots per 

pole per phase” 

IC number of parallel circuits (integer variable) 

IDM multiple of IZY 

IIQQ numerator of fraction ’’slots per pole per phase” 

IPN number of phases (set equal to 3) 

IPX number of poles (integer variable) 

IQQ number of stator slots (integer variable) 

IZY product of number of poles and number of phases 

P number of poles (real variable) 

PBA phase belt angle, deg 

PC number of parallel circuits (real variable) 

PF pitch factor 

PN number of phases (set equal to 3) 
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Subroutine WDGFCT 


QN number of stator slots per pole per phase 

QS number of stator slots (real variable) 


WDGFCT subroutine name 

WDGPCH stator winding pitch expressed as decimal fraction, per unit 
YY slots spanned per armature coil (number slots between coil sides plus 1) 


Subroutine CMBNTN 

Definitions of those variables that are not Listed are the same as in the main pro- 
gram. 

A 3* FLOAT(I)*P, where I = 1,2,3,..., 1000 

CMBNTN subroutine name 
D ABS(QS-F) 

F number of rotor bars 

FF rotor skew, expressed as fraction of rotor circumference, necessary to 

eliminate certain undesirable characteristics in torque-speed curve 

I index 

II index 

K indicator (if K = 1 the slot combination is found to be undesirable; the sub- 

routine will then search for an alternate number of rotor slots) 

L F (L is an integer variable) 

M an indicator showing seriousness of an undesirable slot combination (M = 1 

is most serious; M = 3 is least serious) 

NB number of rotor slots 

P number of poles 

QS number of stator slots 

X constant (1. OE- 15) 
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Main program: INDMTR 


Read one data card 

If card has "M" "Tlf card lacks "M" 


in card column 1 

(a) Read stator 
and rotor 
material data 

(b) Read windage 
loss data 

{c> Read motor 
design data ^ 


in card column 1, 
read motor design 
data 


Calculate various dimensions, 
areas, Carter coefficients, and 
equivalent circuit parameters 



Write first two pages of output 

No-load magnetic calculations 

(i terate on RO an d XO ) 

Write results of magnetic 
calculations (page 3 of output) 

Scale windage loss from 
reference to design conditions 


Write results of windage loss 
scaling (page 3 of output) 

Perform equivalent circuit 
analysis and print results 
(pages 4 and 5 of output) 


Plot torque-speed curve 
(page 6 of output) 


Go back to start of program 


/ 

/ 

✓ 

Called once for stator 
and once for rotor slots 


Subroutine SLOTS 


Calculates stator and 
rotor slot dimensions 
A and slot leakage 
/ I permeance ratios 




/ 


/ 


Subroutine CMBNTN 


Checks that number 
of rotor slots is 
compatible with 
number of stator slots 


Subroutine WDGFCT 


\ Called repeatedly 
\during iteration 

\ '"n 

Called repeatedly 
during iteration 
(SLIP =0) 

\ 

\ 

Called once for \ 
each value of slip 


Calculates distribution! 
and pitch factors for 
stator winding 


Subroutine MAGNET 


Performs magnetic 
circuit calculations 


Subroutine CIRCT 








Does equivalent 
J circuit analysis for 
^ one value of slip 


\ 








\ 


\ 


Subroutine PLOTXY 

Does actual plotting; 
is not part of INDMTR 
program (ref. 3) 


Figure 2. - Simplified flow chart of induction motor computer program. 
















(a) Slot type 1 - rectangular, open. 



(b) Slot type 2 - trapezoidal, open. 



(c) Slot type 3 - rectangular, partially closed. 



(e) Slot type 5 - round. 




partially closed. 


Figure 4. - Allowable rotor and stator slot types with dimensions needed to calculate slot permeance ratio. 
(Symbols shown are those used in subroutine SLOTS. To change the symbols to those used In main 
program, replace each X with S for stator slots or each X with R for rotor slots. For other dimensions 
see figs. 12 and 14) 
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X0 = 0.5*X0AG 

WO = (WSYOKE + WSTOTH)«WFE«3. 0 
R0 = 5.«VW1/W0 
S = 0. 




Calculate total flux and flux per 
pole as follows: 

P-V2 


FTOTAL = 


FPOLE = 


1. 4141*N*F’»KPS«KDS 
0. 637*FTOTAL 


Call CIRCT to compute V2 and IMAG 


> 


Call MAGNET to compute all flux 
densities (BSY and BST in 
particular) and all ampere-turn 
drops and the total ampere-turn 
drop ATTOT 



Compute new values of WO and RO: 


WSYOKE* 







End of no-load 
magnetic calculations 


Figure 5. - Flow chart of no-load magnetic calculations. 


[These statements 
are executed once 
for each data set 


The statements below 
this line are executed 
once for each motor 
design 



Figure 6. - Flow chart of synchronous windage loss calculation. 
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Figure 9. - Material deck for M-19 silicon steel, 




Figure 10. - Magnetization curve for M-19 steei. 



Actual core-loss 

curve 

Core-loss curve 

as reconstructed 
by program 



Figure 11. - Core loss as function of frequency for M-19 steel (0.014-in. -thick lami- 
nations) at 64.5 kilolines per square inch. 
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(c) Slot type 3 - rectangular, partially closed. (d) Slot type 4 - trapezoidal, partially closed. 



(e) Slot type 5 - round. (f) Slot type 6 - trapezoidal with rounded bottom, 

partially closed. 

Figure 12. - Slot dimensions allowable as program input (Symbols shown are those used in subroutine SLOTS. 
To change to symbols used in main program, replace each X with S for stator slots or each X with R for rotor 
slots. The shaded area is CAREA in the notation of subroutine SLOTS. In the main program the shaded area 
Is called S CAREA for stator slots and SB for rotor slots. Where shaded area is shown in two halves, it is 
assumed each half is CAREA/2. For other slot dimensions see figs. 4 and 14. ) 
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(a) Slot type 2 - trapezoidal open. 



(b) Slot type 4 - trapezoidal, partially closed. 



partially closed. AR + AS = CAREA/2. 

Figure 14. - Slot dimensions used in subroutine SLOTS. 
(For other slot dimensions see figs. 4 and 12. ) 


A-Langley, 1976 


E-8486 
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